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ABSTRACT 


HAGERMAN,  LISA  MARIE.  Regional  Analysis  of  Non-methane  Hydrocarbons  and 
Meteorology  of  the  Rural  Southeast  United  States.  (Under  the  direction  of  Viney  P. 
Aneja). 

Measurements  of  non-methane  hydrocarbons,  as  well  as  ozone,  meteorological  and 
trace  gas  data,  were  made  at  four  rural  sites  located  within  the  southeastern  United  States 
as  a  part  of  the  Southern  Oxidants  Study.  Fifty-six  C2-C10  hydrocarbons  were  collected 
from  1200-1300  local  time,  once  every  six  days  from  September  1992  through  October 
1993.  The  measurements  were  made  in  an  effort  to  enhance  the  understanding  of  the 
behavior  and  trends  of  ozone  and  other  photochemical  oxidants  in  this  region.  The  light 
molecular  weight  alkanes  (ethane,  propane,  n-butane,  iso-butane),  ethene  and  acetylene 
display  a  seasonal  variation  with  a  winter  maximum  and  summer  minimum.  Isoprene  was 
virtually  non-existent  during  the  winter  at  all  sites,  and  averaged  from  9.8  ppbC 
(Yorkville,  GA)  to  21.15  ppbC  (Centreville,  AL)  during  the  summer.  The  terpene 
concentration  was  greatest  in  the  summer  with  averages  ranging  between  3.19  ppbC 
(Centreville,  AL)  to  6.38  ppbC  (Oak  Grove,  MS),  but  was  also  emitted  during  the  winter 
months,  with  a  range  of  1.25  to  1.9  ppbC  for  all  sites.  Propylene-equivalent 
concentrations  were  calculated  to  account  for  differences  in  reaction  rates  between  the 
hydroxyl  radical  and  individual  hydrocarbons,  and  to  thereby  estimate  their  relative 
contribution  to  ozone,  especially  in  regards  to  the  highly  reactive  biogenic  compounds 
such  as  isoprene.  It  was  calculated  that  biogenics  represent  at  least  65%  of  the  total  non¬ 
methane  hydrocarbon  sum  at  these  four  sites  during  the  summer  season  when  considering 
propylene-equivalent  concentrations.  An  ozone  episode  which  occurred  from  July  20  to 
July  24  1993  was  used  as  an  example  to  show  ozone  profiles  at  each  of  the  sites,  and  to 
show  the  effect  of  synoptic  meteorology  on  high  ozone  by  examining  NOAA  daily 
weather  maps  and  climatic  data.  Relationships  between  meteorological  variables  such  as 
temperature,  relative  humidity,  solar  radiation,  and  trace  gases  such  as  ozone,  SO2,  CO 
and  NOy  were  also  examined.  A  multiple  regression  on  ozone  at  the  Yorkville  site  found 


temperature,  relative  humidity  and  NOy  to  be  statistically  significant,  with  an  R-square  of 
0.66. 

A  multiple  regression  using  meteorological  and  trace  gas  data  as  input  parameters 
and  ARMA  (autoregressive  moving  average)  time  series  errors  was  applied  to  model  daily 
average  ozone,  using  data  collected  from  June  1  through  August  31  1992  at  Yorkville,  a 
rural  site  located  in  Pauldin  County,  Georgia.  This  statistical  model  was  then  applied  to 
the  same  site  for  the  1993  summer,  and  the  predicted  and  observed  ozone  values  were 
compared  using  statistical  tests.  It  was  found  that,  while  a  good  regression  model  was 
fitted  for  the  1992  ozone  data,  the  same  model  tended  to  underpredict  observed  ozone 
concentrations  for  the  1993  summer.  Ozone  was  found  to  be  statistically  higher  in  1993 
than  in  1992  (with  a  1993  average  of  63.5  ±  16.6  ppb  and  a  1992  average  of  50.25  ±  14.8 
ppb),  and  this  appears  to  be  explained  by  the  synoptic  meteorology  characterizing  the 
1993  summer,  which  was  conducive  to  high  ozone  formation.  A  multiple  regression 
model  with  time  series  errors  was  not  found  to  be  adequate  in  predicting  ozone 
concentrations,  because  other  factors  affecting  ozone  levels,  such  as  synoptic  meteorology 
and  atmospheric  chemistry,  also  need  to  be  considered. 
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INTRODUCTION 


It  has  long  been  established  that  non-methane  hydrocarbons  play  an  important  role 
as  contributors  to  ozone  and  other  secondary  photochemical  pollutants.  In  the  1950s, 
Haagen-Smit  and  his  co-workers  wrote  a  classic  series  of  papers  which  explained  the 
smog  problem  in  Los  Angeles,  establishing  that  a  major  component  of  photochemical 
smog  was  ozone  formed  from  reactions  between  volatile  organic  compounds  (VOCs)  and 
oxides  of  nitrogen  (NO*)  in  the  presence  of  sunlight.  Since  the  passage  of  the  1970  Clean 
Air  Act  amendments,  regulatory  efforts  to  comply  with  the  National  Ambient  Air  Quality 
standard  for  ozone  have  largely  failed  [NRC  1991].  Ozone  exceedences  continue  to  be  a 
major  problem,  especially  in  the  southeast  region  of  the  United  States. 

Volatile  organic  compounds  (VOCs)  are  organic  compounds  that  exist  in  the 
vapor  phase  at  standard  temperature  and  pressure.  They  are  significant  because  they  react 
with  oxides  of  nitrogen  to  produce  ozone,  which  is  a  principle  constituent  of 
photochemical  smog.  Hydrocarbons  are  a  subset  of  VOCs  that  contain  only  carbon  and 
hydrogen  atoms.  In  atmospheric  chemistry,  hydrocarbons  are  usually  referred  to  as  non¬ 
methane  hydrocarbons  because  methane  is  relatively  unreactive  in  the  atmosphere  (its 
lifetime  with  the  OH  radical  is  approximately  12  years,  assuming  a  concentration  of  7.7  x 
10*  molecule/cm*)[Vaghjiani  and  Ravishankara  1991],  and  is  thought  to  be  not  important 
in  an  air  quality  perspective. 

There  are  four  classifications  of  hydrocarbons.  Alkanes,  also  called  paraffins,  are 
hydrocarbons  which  contain  only  single  bonds  and,  as  a  group,  are  the  generally  the  least 
reactive  with  the  OH  radical.  Alkenes  are  hydrocarbons  with  double  bonds  and  are  also 
called  olefins.  Alkynes  are  hydrocarbons  with  triple  bonds;  of  the  60  sampled 
hydrocarbons  in  this  work,  acetylene  is  the  only  alkyne.  These  three  groups  are 
collectively  called  “aliphatics.”  A  fourth  classification  include  the  aromatics,  or  arenes, 
which  have  a  structural  unit  based  on  the  benzene  ring. 

Hydrocarbons  are  precursors  to  tropospheric  ozone  formation.  Ozone  exists  in 
two  different  regions  of  the  atmosphere.  In  the  stratosphere  at  an  altitude  of  about  10 
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km,  high  ozone  concentrations  absorb  ultraviolet  radiation  of  200-300  nm,  preventing  it 
from  reaching  the  earth’s  surface;  for  this  reason,  it  is  called  “good  ozone”.  In  the 
troposphere,  ozone  is  a  pollution  problem  because  it  is  a  major  component  of 
photochemical  smog.  Excessive  ozone  concentrations  can  cause  eye  and  bronchial 
irritation,  respiratory  disease,  and  damage  to  forests  and  agriculture  [NRC  1991,  3 1]. 

In  the  troposphere,  ozone  is  produced  by  the  photodisassociation  of  NO2,  as 
shown  in  the  three  step  mechanism  below.  The  chemical  equations  involving  the 
equilibrium  of  ozone  and  reactions  between  NO2,  O2,  O3  and  NO  are  called  the 
Photostationary  State; 


NO2  +  hv 

NO  +  0(^P) 

(1) 

0(^P)  +  O2 

-» 

O3 

(2) 

O3+NO 

-> 

NO2  +  O2 

(3) 

Nitrogen  dioxide  (NO2)  is  photodisassociated  into  nitric  oxide  (NO)  and  an  excited  state 
of  oxygen  (0(^P)).  The  excited  oxygen  reacts  with  a  diatomic  oxygen  molecule, 
producing  ozone.  However,  ozone  reacts  with  nitric  oxide,  forming  NO2  and  O2,  closing 
the  cycle  and  resulting  in  no  net  ozone  accumulation. 

In  the  presence  of  volatile  organic  compounds,  RO2  radicals  are  produced  by  the 
reaction  of  the  VOCs  with  OH  radicals  that  also  exist  in  the  ambient  atmosphere.  This 
competition  between  RO2  radicals  and  O3  occurs  for  the  oxidation  of  NO  to  NO2.  The 
following  set  of  chemical  equations  is  a  simple  example  of  a  straight-chain  alkane 
oxidation: 


R1CH2R2  +  OH 

-> 

RiC*HR2  +H2O 

(4) 

RiR2C*H  +  O2 

R1R2CHOO* 

(5) 

R1R2CHOO*  +NO 

RiR2CHO’  +NO2 

(6) 

RiR2CHO’ 

-> 

R2CHO  +  Ri 

(7) 

RiR2CHO’  +  O2 

-> 

R1COR2  +  HO2 

(8) 

Ri  O2 

-> 

RiOO* 

(9) 

RiOO*  +  NO 

-> 

NO2  +  RiO 

(10) 
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RiO’  +  O2 

Ri’CHO  +  HO2 

(11) 

HO2  +  NO 

NO2  +  OH 

(12) 

Net: 

R1CH2R2  +  3N0  +  02^  3NO2  +  R1COR2  +  R’lCHO  +  H2O  (13) 

NO2  +  hv  NO  +  O3  (14) 

In  this  example  Ri  and  R2  stand  for  alkyl  groups,  Ri’  represents  an  alkyl  group  with  one 
carbon  atom  less  than  Ri  (i.e.  RiO*  =  Ri’CH20).  A  free  radical  is  a  species  that  contain 
unpaired  electrons  (i.e.  ROO*  is  a  peroxy  radical  [Wameck  1988,  190]).  Three  major 
things  occur  in  this  reaction  sequence.  First,  a  hydrogen  atom  is  abstracted  from  the 
alkane  by  the  hydroxyl  radical  (OH)  to  form  water  and  an  organic  radical  (equation  4), 
which  combines  with  O2  in  equation  (5)  to  form  the  peroxy  radical.  This  peroxy  radical 
reacts  with  nitric  oxide  to  form  NO2  (eq.6),  and  NO2  is  photolyzed  to  produce  ozone. 

The  by  products  of  this  oxidation  include  a  ketone  (R1COR2)  and  an  aldehyde  (RiCHO), 
together  called  “carbonyls.” 

The  contribution  of  naturally  emitted  volatile  organic  compounds  to  ozone 
formation  has  become  of  greater  concern  within  the  last  decade  [Lamb  et  al.,  1987; 
Trainer  et  al.,  1987;  Chameides  et  al.,  1988].  Isoprene,  a-pinene,  P-pinene  and  limonene 
are  some  examples  of  hydrocarbons  emitted  by  vegetation.  Isoprene  (2-methyl- 1,3- 
butadiene)  is  primarily  emitted  from  deciduous  trees,  while  alpha  and  beta  pinene  are 
emitted  from  coniferous  species,  which  emit  three  times  less  than  deciduous  isoprene 
emitters,  on  a  unit  mass  basis  [Lamb  et  al.,  1987].  Almost  half  of  the  isoprene  emitted  in 
the  United  States  is  from  the  Southeast  region  (48%),  while  the  Northwest  accounts  for 
only  3%  of  emissions  [Lamb  et  al.,  1987].  Isoprene  emissions  are  greatest  during  the 
summer  when  temperatures  are  highest,  and  minimal  in  the  wintertime.  Studies  by  Tingey 
[1981]  conducted  on  live  oak  and  slash  pine  using  controlled  environmental  growth 
chambers  found  that  temperature  increased  both  isoprene  and  monoterpene  emissions. 
Tingey  also  found  that  isoprene  was  emitted  by  the  application  of  simulated  daylight,  and 
the  emission  rate  increased  with  light  intensity.  Monoterpene  emissions  were  not  directly 
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affected  by  light.  Natural  emissions  of  isoprene,  the  monoterpenes,  and  other  biogenic 
VOCs  are  a  concern  because  of  the  high  reaction  rate  between  these  compounds  with 
ozone  and  the  OH  radical,  and  therefore,  the  strong  potential  to  act  as  precursors  to 
tropospheric  ozone  formation.  Measurements  suggest  that  isoprene  and  other  biogenics 
result  in  high  ozone  concentrations  in  urban  areas  affected  by  high  NOx  concentrations,  as 
well  as  contributing  to  relatively  high  ozone  in  rural  areas  [NRC  1991,  8]. 

Summertime  ozone  in  the  South  is,  on  average,  among  the  highest  in  the  United 
States  [SOS  1995,  1 1],  One  reason  for  this  is  due  to  the  stagnant  conditions  during  the 
summer,  which  inhibit  dispersion  of  pollutants  and  allows  a  steady  buildup  of  pollutants 
and  ozone  precursors.  These  stagnating  high  pressure  systems  often  arise  from  the  stalling 
of  a  continental  high-pressure  system  over  the  Appalachian  Mountains  [SOS  1993,  12]. 
Another  characteristic  which  is  conducive  to  high  ozone  in  this  region  is  the  dense 
vegetation,  which  results  in  large  emissions  of  isoprene  and  other  natural  hydrocarbons 
[Lamb  et  al.,  1987].  These  natural  hydrocarbons  are  highly  reactive,  and  model  studies 
indicate  that  these  natural  emissions  can  significantly  affect  both  urban  and  rural  ozone 
levels  [Chameides  et  al.,  1988;  Trainer  et  al.,  1987].  There  have  been  various  studies  done 
involving  the  measurement  of  background  levels  of  hydrocarbons  in  the  atmosphere  of 
rural  continental  sites  [Jobson  et  al.,1994;  Hov  et  al.,  1991;  Colbeck  and  Harrisson,  1985; 
Greenberg  and  Zimmerman,  1984],  however,  prior  to  this  study,  an  analysis  of  rural 
hydrocarbons  in  the  southeast  United  States  has  not  been  done  on  a  regional  scale.  In  this 
study,  fifty-six  non-methane  hydrocarbons  were  analyzed  from  four  rural  sites  located  in 
the  southeastern  United  States,  along  with  ozone,  various  trace  gases,  and  meteorological 
variables.  Propylene-equivalent  concentrations  were  calculated  to  evaluate  the 
contribution  of  hydrocarbons  to  ozone,  especially  the  highly  reactive  biogenic  compounds. 
Relationships  between  ozone,  isoprene  and  various  meteorological  parameters  were  also 
examined. 
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SECTION  I 


REGIONAL  CHARACTERIZATION  OF  NON-METHANE 
HYDROCARBONS  IN  THE 
RURAL  SOUTHEAST  UNITED  STATES 


Abstract 

Concentrations  of  non-methane  hydrocarbons,  as  well  as  ozone,  meteorological 
and  trace  gas  data,  were  measured  at  four  rural  sites  located  within  the  southeastern 
United  States  as  a  part  of  the  Southern  Oxidants  Study.  Fifty-six  C2-C10  hydrocarbons 
were  collected  from  1200-1300  local  time,  once  every  six  days  from  September  1992 
through  October  1993.  The  measurements  were  made  in  an  effort  to  enhance  the 
understanding  of  the  behavior  and  trends  of  ozone  and  other  photochemical  oxidants  in 
this  region.  The  light  molecular  weight  alkanes  (ethane,  propane,  n-butane,  iso-butane), 
ethene  and  acetylene  display  a  seasonal  variation  of  a  winter  maximum  and  summer 
minimum.  Isoprene  was  virtually  non-existent  during  the  winter  at  all  sites,  and  averaged 
from  9.8  ppbC  (Yorkville,  GA)  to  21.15  ppbC  (Centreville,  AL)  during  the  summer.  The 
terpene  concentration  was  greatest  in  the  summer  with  averages  ranging  3.19  ppbC 
(Centreville,  AL)  to  6.38  ppbC  (Oak  Grove,  MS),  but  was  also  emitted  during  the  winter 
months,  v^th  a  range  of  1.25  to  1.9  ppbC  for  all  sites.  Propylene-equivalent 
concentrations  were  calculated  to  account  for  differences  in  reaction  rates  between  the 
hydroxyl  radical  and  individual  hydrocarbons,  and  to  thereby  estimate  their  relative 
contribution  to  ozone,  especially  in  regards  to  the  highly  reactive  biogenic  compounds 
such  as  isoprene.  It  was  calculated  that  biogenics  represent  at  least  65%  of  the  total  non¬ 
methane  hydrocarbon  sum  at  these  four  sites  during  the  summer  season  when  considering 
propylene-equivalent  concentrations.  An  ozone  episode  which  occurred  from  July  20  to 
July  24  1993  was  used  as  an  example  to  show  ozone  profiles  at  each  of  the  sites,  and  to 
show  the  effect  of  synoptic  meteorology  on  high  ozone  by  examining  NOAA  daily 
weather  maps  and  climatic  data.  Relationships  between  meteorological  variables  such  as 
temperature,  relative  humidity,  solar  radiation,  and  trace  gases  such  as  ozone,  SO2,  CO 
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and  NOy  were  also  examined.  A  multiple  regression  on  ozone  at  the  Yorkville  site  found 
temperature,  relative  humidity  and  NOy  to  be  statistically  significant,  with  an  R-square  of 
0.66. 

1.  Introduction 

It  has  long  been  established  that  non-methane  hydrocarbons  play  an  important  role 
as  precursors  to  ozone  and  other  secondary  photochemical  pollutants.  Ozone  is  formed 
from  reactions  between  volatile  organic  compounds  (VOCs)  and  oxides  of  nitrogen  (NO*) 
in  the  presence  of  sunlight.  Since  the  passage  of  the  1970  Clean  Air  Act  amendments, 
regulatory  efforts  to  comply  with  the  National  Ambient  Air  Quality  standard  for  ozone 
have  been  inadequate  [NRC  1991,  4;  Dimitriades  1989].  Ozone  exceedences  continue  to 
be  a  major  problem,  especially  in  the  southeast  region  of  the  United  States.  Studies  have 
shown  that  the  Southeast  is  a  region  where  high  concentrations  of  ozone  accumulate  in 
both  rural  and  urban  areas  [SOS  1995,  iv].  The  contribution  of  naturally  emitted  volatile 
organic  compounds  to  ozone  formation  in  both  urban  and  rural  areas  has  become  of 
greater  concern  within  the  last  decade  [Lamb  et  al.,  1987].  Measurements  of  biogenically 
emitted  VOCs  such  as  isoprene  suggest  that  these  compounds  contribute  to  high  ozone 
concentrations  in  urban  areas  affected  by  NO*  [Trainer  et  al.,  1987;  Chameides  et  al., 
1988;  NRC  1991,  8].  Various  studies  have  involved  measuring  isoprene  and  other 
ambient  hydrocarbon  concentrations  in  rural  or  remote  sites  [Andronache  et  al.  1994; 
Chameides  et  al.  1992;  Colbeck  and  Harrisson,  1985;  Greenberg  and  Zimmerman,  1984; 
Rasmussen  and  Khalil,  1988;  Sexton  and  Westberg,  1984].  Other  studies  have  reported 
the  seasonal  variations  of  hydrocarbons  in  continental  air  [Boudries  et  al.,  1994;  Jobson  et 
al.,  1994;  Hov  et  al.,  1991;  Rudolph  et  al,  1989;  Tille  et  al,  1985]. 

In  this  study  we  (1)  compare  C2-C10  hydrocarbons  during  maximum  photochemical 
activity  on  a  regional  scale,  (2)  analyze  the  contribution  of  rural  hydrocarbons  using 
propylene-equivalent  concentrations,  especially  in  relation  to  isoprene,  and  (3)  examine 
the  relationship  between  ozone  and  reactive  nitrogen  (NOy) ,  and  between  ozone  and 
meteorological  variables  such  as  temperature,  relative  humidity,  solar  radiation,  and 
ozone. 
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2.  Experiment 
2.1.  Site  Description 

The  hydrocarbon,  trace  gas  and  meteorological  data  was  collected  from  four  rural 
sites  within  the  SOS-SCION  network  (Southern  Oxidants  Study-Southeastern  Consortium 
Intermediate  Oxidant  Network)  located  in  the  southeast  United  States.  The  SCION 
network  was  established  to  describe  how  ozone  precursor  concentrations  vary  during  the 
year  in  different  regions  of  the  Southeast.  The  location  of  the  sites  used  for  this  paper 
include  Centreville,  Alabama;  Oak  Grove,  Mississippi;  Yorkville,  Georgia;  and  Candor, 
North  Carolina.  Figure  1  illustrates  the  location  of  the  sites. 

The  Centreville  site  is  located  in  Bibb  County,  Alabama  (32°90’N,  87°23’W),  in  a 
rural  area  representative  of  the  transitional  nature  of  the  region  between  the  lower  coastal 
plain  and  Appalachian  highlands,  at  an  elevation  of  136  m  mean  sea  level  (msl).  Sources 
of  anthropogenic  emissions  located  vdthin  a  1 10  km  radius  of  the  sampling  site  include  the 
cities  of  Montgomery,  Birmingham,  and  Tuscaloosa.  This  site  is  located  in  a  large  field 
approximately  180  m  from  a  NOAA  (National  Oceanic  and  Atmospheric  Administration) 
weather  radar  station. 

The  Oak  Grove  site  (30°99’N,  88°93’W)  is  located  in  the  Desoto  National  Forest 
in  Perry  County,  Mississippi  at  an  elevation  85  m  msl.  This  site  is  located  in  a  rural  area 
representative  of  the  lower  coastal  plain.  It  is  moderately  forested  with  a  canopy  at 
approximately  12  m,  and  predominantly  consists  of  conifers.  The  forested  areas  are 
interspersed  with  cultivated  farm  land.  The  site  is  located  in  a  large  field  approximately  46 
m  from  the  nearest  row  of  brush  and  152  m  from  the  intersection  of  two  dirt  roads.  The 
area  immediately  surrounding  the  site  is  mowed  on  a  regular  basis  and  is  surrounded  by 
cultivated  farm  land.  The  nearest  residence  visible  from  the  site  is  located  approximately  a 
quarter  of  one  mile  away.  It  is  situated  off  of  state  route  29  and  approximately  40  km 
southeast  of  Hattiesburg,  Mississippi. 
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Figure  1.  Map  of  sampling  sites 
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The  Yorkville  site  (33°55’4r’N,  85°02’46”W)  is  located  in  Pauldin  County, 
Georgia,  at  an  elevation  approximately  400  m  above  sea  level.  The  site  is  situated  in  a 
rural  area  representative  of  the  southern  highlands,  consisting  of  hardwood  forests 
interspersed  with  open  pasture  and  tilled  farmland.  The  site  is  off  route  278  and  is 
approximately  72  km  west  of  Atlanta,  and  48  km  west  of  a  power  generating  station. 

The  Candor  site  (35.26°N,  79.84°W,  197  m  msl  elevation)  is  located  in  the  Central 
Piedmont  region  of  North  Carolina  on  the  eastern  border  of  the  Uwharrie  National  Forest. 
The  sampling  site  is  located  in  an  open  field  approximately  1200  m^,  and  the  field  is 
surrounded  by  forests  mixed  with  deciduous  and  coniferous  trees.  Sources  of 
anthropogenic  pollution  located  within  a  120  km  radius  of  the  sampling  site  include  the 
urban  areas  ofRaleigh-Durham,  Greensboro,  Winston-Salem,  and  the  junction  between  I- 
40  and  1-85,  which  are  all  situated  to  the  north  and  northeast  of  the  site. 

2.2.  Data  Collection  and  Analysis 

The  hydrocarbon  data  used  in  this  paper  was  sampled  from  September  1992 
through  October  1993,  and  includes  data  sampled  during  June  1992  at  the  Candor  site. 

The  Centreville  and  Candor  sites  each  include  49  total  observations;  Oak  Grove  and 
Yorkville  include  61  and  55  observations,  respectively.  The  sampling  days  used  for  the 
seasonal  averages  at  each  site  are  listed  at  the  end  of  Tables  la-d.  Time  integrated  air 
samples  were  collected  in  6  liter  SUMMA  electropolished  stainless  steel  canisters  from 
1200  to  1300  local  time,  once  every  six  days.  This  sampling  frequency  was  chosen  so  that 
each  day  of  the  week  would  be  represented  in  the  study.  The  hydrocarbon  samples  were 
collected  in  evacuated  canisters  by  opening  the  canister  and  allowing  the  internal  pressure 
to  reach  ambient  pressure.  C2-C10  hydrocarbons  were  analyzed  at  the  University  of  Miami 
using  a  Hewlett  Packard  HP  5890II  gas  chromatograph  equipped  with  a  cryogenic  cooling 
option  and  flame  ionization  detection.  Data  reduction  was  accomplished  using  HP  3365 
Chemstation  II  software  on  PC-DOS  based  personal  computers.  The  automatic  air 
concentrator  used  was  a  modified  Entech  2000  (Entech  Laboratory  Automation,  Simi 
Valley,  CA).  The  detection  limit  was  0. 1  ppbC  (parts  per  billion  carbon)  with  a 
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reproducibility  of  30%.  A  detailed  description  of  the  GC  analysis  is  published  elsewhere 
[Farmer  et  al.,  1994], 

Measurements  of  NO,  NOy  (the  sum  of  the  reactive  odd  nitrogen  species,  NOy  = 
NO  +  NO2  +  organic  nitrates  +  inorganic  nitrates),  SO2,  CO,  and  O3  as  well  as 
meteorological  parameters  such  as  temperature,  relative  humidity,  solar  radiation, 
barometric  pressure,  and  wind  speed  and  direction  were  made  every  day  at  fifteen  minute 
intervals  throughout  the  year.  NO  and  NOy  were  measured  with  the  TECO  42S  (Thermo 
Environmental  Instruments  Inc.)  chemiluminescent  high  sensitivity  analyzer.  For  the 
Centreville,  Oak  Grove  and  Yorkville  sites,  air  samples  for  the  continuous  gas  monitoring 
equipment,  excluding  the  non-methane  hydrocarbons,  were  collected  through  0.25”  teflon 
tubing,  each  instrument  equipped  with  a  dedicated  teflon  line  and  particulate  filter  located 
at  the  intake.  The  NOy  converter  for  the  TECO  42S  was  located  within  the  intake  dome 
and  operated  at  350°C.  Cylinders  containing  gas  standards  for  NO,  NOy,  SO2  and  CO 
were  present  along  with  a  TECO  146  dynamic  gas  calibrator  to  provide  for  calibration  and 
zero  and  span  checks.  The  site  was  equipped  with  a  zero  air  generating  system.  Zero  air 
for  the  O3,  NO/NOy  and  SO2  monitors  was  generated  by  passing  ambient  air  through  a 
series  of  canisters  containing  purafill,  activated  charcoal  and  brominated  charcoal.  Zero 
air  for  the  CO  monitor  was  generated  by  passing  ambient  air  through  a  palladium 
converter.  The  trace  gas  (except  ozone)  and  meteorological  data  used  for  this  work  was 
an  average  of  the  1200-1300  data  from  the  corresponding  hydrocarbon  sampling  days. 

The  daily  maximum  value  was  used  for  ozone,  which  generally  occurred  mid-afternoon 
(-1500).  The  data  used  for  the  multiple  regression  analysis,  including  CO,  SO2,  NOy  and 
ozone,  plus  the  meteorological  variables,  consisted  of  daily  averages  of  the  1000-1600 
period  from  June  1  1993  to  August  31  1993  (92  observations). 
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3.  Results  and  Discussion 


3.1.  Measurements  of  C2-C10  Speciated  Hydrocarbons 

Tables  la-d  summarize  the  seasonal  average,  median,  standard  deviation  and  range 

of  the  fifty-six  C2-C10  compounds  sampled  at  each  site.  Isobutene  and  1 -butene  were 
added  together  because  of  problems  with  coelution,  as  well  as  m-xylene  and  p-xylene.  For 
the  most  part,  the  autumn  through  summer  seasons  represent  the  data  in  chronological 
order:  autumn  is  associated  with  September  through  November  1992,  winter  includes 
December  1992  through  February  1993,  spring  includes  March  through  May  1993,  and 
summer  includes  June  through  August  1993.  However,  the  autumn  category  also  includes 
sampling  days  from  September  and  October  1993,  and  the  summer  category  includes  three 
days  from  June  1992  at  the  Candor  site.  The  sampling  days  used  for  the  seasonal  averages 
are  listed  at  the  end  of  Tables  la-d. 

During  the  wintertime,  all  sites  had  the  same  top  four  dominant  compounds  in  the 
following  order;  propane,  n-butane,  ethane,  and  isopentane,  with  the  exception  of  the 
Yorkville  site,  which  had  isopentane  and  ethane  switched  around  in  ranking.  The 
compounds  following  the  ones  listed  above  were  within  the  top  10  most  abundant 
species:  acetylene,  n-pentane,  ethene,  isobutane,  and  benzene.  Toluene  was  also  within 
the  top  10  at  all  sites  except  at  Oak  Girove,  where  it  ranked  1 1th.  The  individual  C2-C5 
alkanes  (except  cyclopentane)  dominated  the  list  of  most  abundant  compounds  for  all  four 
sites  during  the  winter.  During  the  summer,  the  10  most  abundant  compounds  were 
highly  variable  among  the  sites.  Isoprene  was  the  dominant  hydrocarbon  during  the 
summer  at  all  but  the  Oak  Grove  site,  where  n-pentane  was  the  most  abundant 
hydrocarbon  (16.51  ±  20.97  ppbC,  median  =  7.68  ppbC)  with  isoprene  immediately 
following  with  a  concentration  of  1 1  ±  4.2  ppb.  The  compounds  isoprene,  propane, 
isopentane,  2-methylpentane,  and  styrene  consistently  appeared  in  the  top  10  most 
abundant  hydrocarbons  among  all  four  sites,  though  not  necessarily  in  that  order. 

Figures  2-4  shows  the  seasonal  averages  of  various  compounds  and  compound 
sums  for  each  site.  In  Figure  2,  ethane,  propane,  n-butane,  isobutane,  ethene  and 
acetylene  display  a  distinct  seasonal  variation  with  maximums  occurring  during  the  winter. 
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Table  la.  Seasonal  averages  of  hydrocarbon  compound  concentrations  in  ppbC  at  Centreville,  Alabama _ 
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Table  Ic.  Seasonal  averages  of  hydrocarbon  compound  concentrations  in  ppbC  at  Yorkville,  Georgia _ 

Autumn  _ Winter  _ Spring _  Summer 

Compound  Mean  Medn  StD  Max  Min  Mean  Medn  StD  Max  Min  Mean  Medn  StD  Max  Min  Mean  Medn  StD  K 
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Table  Id.  Seasonal  averages  of  hydrocarbon  compound  concentrations  in  ppbC  at  Candor,  North  Carolina 


Methylcyclohexane 


Table  Id.  Seasonal  averages  of  hydrocarbon  compound  concentrations  in  ppbC  at  Candor,  North  Carolina _ 

Autumn _  _ Winter _  _  Spring  Summer 

_ Compound _ Mean  Medn  StD  Max  Min  Mean  Medn  StD  Max  Min  Mean  Medn  StD  Max  Min  Mean  Medn  StD 
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This  seasonal  variation  of  the  paraflSns  and  acetylene  is  consistent  with  the  literature.  A 
study  by  Jobson  et  al,  [1994]  found  that  at  a  remote  boreal  site  in  Canada,  alkane 
compounds  and  acetylene  concentrations  displayed  a  winter  maximum  and  summer 
minimum;  Penkett  et  al.,  [1993]  observed  seasonal  trends  in  hydrocarbon  concentrations 
in  air  over  the  North  Atlantic  Ocean.  Seasonal  variations  of  atmospheric  hydrocarbons 
were  also  measured  in  western  France  [Boudries  et  al.  1994]  and  at  a  rural  site  in  Norway 
[Hov  et  al.  1991].  The  winter  maximum  and  summer  minimum  of  the  lower  alkanes  and 
acetylene  has  been  attributed  to  hydroxyl  chemistry  and  the  seasonal  abundance  of  the  OH 
radical  [Jobson  et  al.l994,  Penkett  et  al.  1993,  Boudries  et  al.  1994,  Lightman  et  al.  1990, 
Spivakovsky  et  al.  1990].  The  seasonal  variation  in  hydrocarbon  source  strengths,  and 
differences  in  atmospheric  behavior  such  as  increased  convection  and  vertical  mixing  in  the 
summer,  or  differences  in  air  mass  climatology  with  season,  also  play  a  role  in  the 
hydrocarbon  seasonal  variation  [Jobson  et.al  1994]. 

The  sum  of  C2-C10  hydrocarbons  (Figure  3)  did  not  vary  greatly  between  the 
autumn,  winter  and  spring,  ranging  from  approximately  30  to  50  ppbC  for  these  three 
seasons.  However,  the  summer  period  shows  much  higher  concentrations,  accounted  for 
primarily  by  increased  concentrations  of  the  biogenic  sums. 

The  paraffin  sum  displays  a  seasonal  pattern  with  a  winter  maximum  for  all  but  the 
Oak  Grove  site,  which  displays  its  maximum  during  the  summer.  This  can  be  attributed  to 
n-pentane  (16.51  ±  20.97  ppbC,  median  =  7.68  ppbC),  which  is  almost  8  times  greater 
than  the  second  highest  value,  2. 15  ±  2.28  ppbC  at  the  Yorkville  site,  as  shown  in  Figure 
4.  During  the  summer.  Oak  Grove  had  some  unusually  large  values  of  n-pentane.  Out  of 
12  observations,  2  days  had  values  of-  57  ppbC,  one  day  measured  29  ppbC  and  two 
days  measured  approximately  15  ppbC.  The  reason  for  these  high  n-pentane  values  are 
uncertain;  however,  auto  emissions  is  ruled  out  as  a  possible  source  because  of  the  low 
acetylene  values. 

The  olefins  sum,  which  does  not  include  the  biogenic  hydrocarbons  isoprene,  a- 
pinene,  P-pinene  and  limonene,  shows  little  variation  throughout  the  year,  ranging 
between  ~  4  and  7  ppbC  among  the  sites.  During  the  summer,  however,  olefin 
concentrations  at  Yorkville  are  higher  than  the  other  three  sites  (10.72  ppbC).  Studies 
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Autumn  Winter  Spring  Summer 


Autumn  Winter  Spring  Summer 


Autumn  Winter  Spring  Summer 


Figure  2  continued.  Seasonal  averages  of  select  C2-C4  hydrocarbons 


have  shown  that  there  is  no  distinct  seasonal  trend  for  olefins  at  continental  remote  sites 
unaffected  by  anthropogenic  sources  [Jobson  et  al,  1994],  and  very  little  seasonal 
variation  in  the  free  troposphere  over  the  Atlantic  [Penkett  et  al.,  1993].  Hov  et  al  [1991] 
found  a  seasonal  trend  for  ethene  and  propene  with  a  late  January  maxima  and  a  secondary 
maxima  during  July- August,  but  attributed  this  summer  maxima  to  local  release  from  areas 
upwind  of  the  site.  Yorkville  is  affected  by  two  potential  emission  sources;  the  city  of 
Atlanta,  located  approximately  45  miles  (72  km)  southeast  of  the  Yorkville  site,  and  a 
large  power  generating  station  located  approximately  30  miles  (48  km)  east  of  Yorkville 
[Kirk  1996].  The  summer  maxima  for  the  olefins  at  Yorkville  may  be  due  to  transport 
fi-om  these  anthropogenic  emission  sources. 

While  Oak  Grove  and  Yorkville  display  summer  aromatic  concentrations 
approximately  two  times  higher  than  during  the  other  three  seasons,  the  Candor  site 
displays  extraordinarily  large  concentrations  of  the  aromatic  sums  during  both  the  autumn 
and  summer  seasons.  The  difference  between  the  high  aromatic  concentration  during 
autumn  at  Candor  (19.75  ppbC)  and  the  other  three  sites  (~  5-8  ppbC)  is  accounted  for  by 
benzene  (12.77  ±  8.35  ppbC).  The  difference  between  the  high  summer  aromatic  sum  at 
Candor  and  the  Yorkville  site  is  due  primarily  to  1,2,4-trimethylbenzene  (6.79  ±  4.66), 
along  with  the  compounds  styrene  (2.75  ±  1.28)  and  toluene  (2.94  ±  1.97).  Figure  4 
shows  the  compounds  n-pentane,  benzene,  and  1,2,4-trimethylbenzene,  which  exhibit 
particularly  high  concentrations.  The  reason  for  the  high  levels  of  these  particular 
compounds  is  unclear.  Surface  wind  direction  appears  to  have  a  random  relationship  with 
excessively  high  concentrations  of  n-pentane,  benzene,  and  1,2,4-trimethylbenzene. 

Table  2  summarizes  selected  hydrocarbons  measured  during  the  summer  at  various 
rural  and  remote  sites.  It  can  be  seen  that  hydrocarbon  concentrations  at  the  Centreville 
site  are  comparable  to  the  rural  Norway  and  Maine  sites  listed  in  Table  2.  The  Fraserdale 
site  in  Canada  appears  to  be  the  least  affected  by  anthropogenic  sources,  given  its  lower 
concentrations  of  acetylene  and  the  other  lower  molecular  weight  hydrocarbons. 


24 


ppbC 


n-Pentane 


Autumn 


Winter 


Spring 


Benzene 


Summer 


□  Centrevilie 

□  Oak  Grove 
■  Yorkville 
B  Candor 


□  Centrevilie 

□  Oak  Grove 

■  Yorkville 

■  Candor 


Autumn 

1  1 

Winter  Spring 

Summer 

ppbC 

T 

1 ,2,4-Trimethylbenzene 

□  Centrevilie 
□Oak  Grove 
■Yorkville 

□  Candor 


Figure  4.  Hydrocarbons  with  abnormally  high  values 


Table  2.  Average  hydrocarbon  concentrations  (ppbC)  at  various  rural  sites  during  the 
summer  months 


Compound 

Fiaserdale 

Canada’ 

Birkenes 

Norway^ 

Belfast 

Maine^ 

NW^ 

England 

Centreville 

Alabama^ 

Raleigh 

NC® 

Brazif 

Ethene 

1.67 

2.0 

1.6 

0.61 

3.78 

Acetylene 

0.15 

0.51 

<0.5 

0.6 

0.36 

Ethane 

1.64 

3.05 

3.5* 

14.9 

1.74 

4.18 

Propene 

0.82 

0.5 

4.1 

0.73 

0.54 

.93 

Propane 

0.23 

2.01 

2.0 

10.3 

2.60 

9.59 

1.35 

Isobutane 

0.028 

0.70 

0.5 

0.8 

0.77 

0.61 

n-Butane 

0.06 

1.67 

2.0 

1.1 

1.48 

2.04 

0.96 

Isopentane 

0.04 

1.00 

1.0 

2.22 

4.42 

n-Pentane 

0.065 

0.62 

1.0 

5.1 

1.04 

1.97 

<DL 

2-Methylpentane 

<0.5 

3.34 

1.28 

3-Methylpentane 

<0.5 

0.82 

0.70 

n-Hexane 

5.1 

0.68 

0.87 

<DL 

Cis-3-hexene 

<0.5 

0.70 

Benzene 

0.61 

1.32 

3 

Toluene 

1.20 

8.96 

0.84 

*  Geometric  mean 


Reference _ Sampling  Period 


1  Jobson  et  al.,  1994 

July  to  September  1990,  1991  and  1992,  collected  midmoming 
(0900  to  1200).  Samples  from  June  to  July  1990  were  collected 
throughout  the  day. 

2Hovetal.,  1991 

June-August  1987 

3  Sexton  and  Westberg,  1984 

June-July  1975 

4  Colbeck  and  Harrison,  1985 

May-July  1983 

5  This  study 

Jime-August  1993,  1200-1300  local  time 

6  Lawrimore  et  al.,  1995  (semi-uiban) 

surface,  August  1993, 0500-0800  EDT 

7  Greenberg  and  Zimmerman,  1984. 

surface,  August  and  September  1979  and  1980 

Table  3.  Summer  averages  of  biogenic  hydrocarbons  (ppbC) 

Location 

Isoprene  a-Pinene  b-Pinene  Limonene 

Raleigh,  NC’ 

2.08 

Brazil^ 

12 

2.7 

Niwot  Ridge,  CO^ 

3,15 

1.4 

0.7 

Candor,  North  Carolina'* 

10.01 

2.25 

1.30 

0.69 

Centreville,  Alabama^ 

21.15 

1.60 

1.36 

0.24 

Oak  Grove,  Mississippi* 

11.19 

2.71 

3.04 

0.63 

Yorkville,  Georgia* 

9.8 

0.71 

2.06 

0.46 

1  Lawrimore  et  al.,  1995  (surface,  August  1993, 0500-0800  EDT) 

2  Greenberg  and  Zimmerman,  1984  (surface,  August  and  September  1979  and  1980) 

3  Greenberg  and  Zimmerman,  1984  (surface,  August-September  and  November  1982) 

4  This  work  (June  1992  and  1993, 1200-1300  local  time) 

5  This  work  (June-August  1993  1200-1300  local  time) 


3.2.  Biogenic  Hydrocarbons 

Isoprene  averaged  approximately  2  ppbC  during  the  autumn  and  spring,  and  was 
virtually  nonexistent  during  the  winter  period  for  all  four  sites  (mean  <0.1  ppbC)  (Figure 
5a).  Isoprene  concentrations  were  highest  during  the  summer,  with  Centreville  having 
concentrations  twice  as  high  as  the  other  three  sites.  The  terpenes,  which  include  the 
naturally  emitted  compounds  a-pinene,  |5-pinene  and  limonene,  also  show  a  seasonal 
distribution  with  lowest  concentrations  in  the  winter  and  highest  in  the  summer  (Figure 
5b).  However,  unlike  isoprene,  terpenes  are  emitted  throughout  the  winter,  the  sum 
ranging  between  1.25  and  1.9  ppbC  for  all  sites.  A  study  done  by  Tingey  [1981]  on  live 
oak  found  that  isoprene  was  emitted  only  in  daylight,  and  given  constant  light  conditions, 
the  emission  rate  is  temperature  dependent.  While  isoprene  is  highly  dependent  on 
temperature  and  virtually  negligible  during  the  winter,  the  terpenes  have  a  small  winter 
abundance.  Terpene  emissions  from  slash  pine  do  not  vary  with  light,  but  emission  rates 
are  log-linearly  related  to  temperature  [Tingey  1981].  Table  3  lists  isoprene  and 
monoterpene  data  from  this  work  and  other  literature. 

Figure  6  shows  the  linear  relationship  between  the  logarithm  of  isoprene  (ppbC) 
and  temperature  (°C)  using  data  collected  from  all  four  sites.  The  plot  includes  all  values 
for  which  isoprene  was  greater  or  equal  to  1  ppbC  (April  through  September).  The 
regression  equation  is 

log(isoprene)  =  -0.67944  +  0.056202T 

with  an  R-squared  value  of  0. 53 .  The  units  of  isoprene  are  in  ppbC.  The  slope  of  the  best 
fit  line  (0.056)  is  lower  than  that  found  by  Jobson  et  al.  [1994],  who  found  a  slope  of 
0.071.  The  regression  equation  found  by  Jobson  in  units  of  ppbv  was 

log(isoprene)  =  -1 .40  +  0.071T  (in  units  of  ppbv) 

log(isoprene)  =  -0.70103  +  0.071T  (converted  to  units  of  ppbC) 

Converting  this  equation  for  units  of  ppbC  changes  only  the  intercept  to  a  value  of 
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Figure  5a.  Seasonal  average  of  isoprene 
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Figure  5b.  Seasonal  average  of  terpene  sum  (a-,  P-pinene,  limonene) 
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Figure  6.  Isoprene  concentration  versus  temperature  for  the  months  April  through  September  1992  and  1993. 


-0.70103  and  the  slope  remains  the  same.  Considering  a  temperature  range  of  18°  to 
35°C,  the  predicted  isoprene  concentrations  using  Jobson’s  regression  equation  was  two 
to  three  times  greater  than  the  predicted  concentrations  using  the  regression  equation 
found  in  this  work  (refer  to  Table  4).  Part  of  this  discrepancy  may  be  due  to  the  fact  that 
in  this  work,  four  sites  were  used  for  the  regression,  versus  Jobson’s  one  site,  which 
increases  the  variability  of  isoprene  concentrations.  To  account  for  this  difference,  a 
regression  of  isoprene  at  only  one  site,  Centreville,  resulted  in  a  higher  R-squared  value 
(0.61)  and  a  slope  very  similar  to  Jobson  (0.066),  but  with  a  lower  intercept.  The 
regression  equation  for  Centreville  turned  out  to  be 

log(isoprene)  =  -0.9029  +  0.066T 

Another  reason  which  may  explain  the  difference  in  regression  equation  terms  is  that  the 
relationship  between  isoprene  and  temperature  is  affected  by  the  type  of  surrounding 
vegetation,  as  well  the  atmospheric  concentration  of  OH  and  ozone.  Other  possible 
factors  which  can  influence  measured  isoprene  concentrations  include  the  time  of  day 
when  samples  were  collected,  and  atmospheric  conditions,  such  as  the  height  of  the 
boundary  layer  or  vertical  mixing  and  turbulence  in  the  atmosphere.  Table  5  lists  the 
correlation  coefficients  between  meteorological  parameters  and  biogenically  emitted 
hydrocarbons.  The  data  for  the  correlations  include  sampling  days  from  all  four  seasons. 
Isoprene  is  much  more  dependent  on  temperature  than  it  is  on  solar  radiation.  The 
correlation  coefficient  between  temperature  and  isoprene  ranges  between  0.62  (Candor)  to 
0.76  (Yorkville).  The  correlation  between  isoprene  and  solar  radiation  is  less  than  0.3  for 
all  sites. 
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Table  4.  Comparison  of  predicted  isoprene  values  using  regression  equations  from  Jobson 
[1994],  the  regression  of  all  sites,  and  the  regression  equation  from  only  Centreville. 


Predicted  Isoprene  (ppbC) _  _ Ratios 


T(°C) 

Jobson 

All  sites 

Centreville 

Jobson/all 

Cnvt/all 

18.00 

3.78 

2.15 

1.93 

1.76 

0.90 

20.00 

5.24 

2.78 

2.61 

1.88 

0.94 

21.00 

6.17 

3.17 

3.04 

1.95 

0.96 

22.00 

7.26 

3.61 

3.54 

2.01 

0.98 

23.00 

8.55 

4.10 

4.12 

2.08 

1.00 

24.00 

10.07 

4.67 

4.80 

2.16 

1.03 

25.00 

11.86 

5.32 

5.59 

2.23 

1.05 

26.00 

13.96 

6.05 

6.50 

2.31 

1.07 

27.00 

16.44 

6.89 

7.57 

2.39 

1.10 

28.00 

19.36 

7.84 

8.81 

2.47 

1.12 

29.00 

22.80 

8.92 

10.26 

2.56 

1.15 

30.00 

26.85 

10.15 

11.94 

2.64 

1.18 

31.00 

31.62 

11.56 

13.90 

2.74 

1.20 

32.00 

37.24 

13.15 

16.18 

2.83 

1.23 

33.00 

43.85 

14.97 

18.84 

2.93 

1.26 

34.00 

51.64 

17.04 

21.93 

3.03 

1.29 

35.00 

60.81 

19.39 

25.53 

3.14 

1.32 

Regression  equations  in  units  of  ppbC: 


Jobson  et  al.,  1994: 
This  work,  all  sites: 
This  work,  Centreville: 


log(isoprene) 

log(isoprene) 

log(isoprene) 


0.70103  +  0.071T 
0,67944  +  0.056202T 
0.9029  +  0.066T 
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Table  5.  Correlation  coefficients  between  biogenically  emitted  hydrocarbons  and 
meteorological  variables.  The  values  in  boldface  indicate  correlations  >0.5  or  <-0.5. 


Centreville 

Temp 

Relhum 

pBARO 

Solar 

Max  03 

Isoprene 

a-Pinene  b-Pinene  Limonene 

TEMP 

1.00 

RELHUM 

0.03 

1.00 

pBARO 

-0.25 

-0.49 

1.00 

SOLRAD 

0.32 

-0.65 

0.21 

1.00 

OzoneMax 

0.51 

-0.35 

-0.07 

0.48 

1.00 

Isoprene 

0.70 

0.02 

-0.07 

0.16 

0.11 

1.00 

a-Pinene 

0.55 

0.32 

-0.23 

-0.13 

-0.05 

0.76 

1.00 

b-Pinene 

0.34 

0.14 

-0.35 

0.01 

0.10 

0.33 

0.39  1.00 

Limonene 

-0.32 

0.26 

-0.10 

-0.12 

-0.16 

-0.27 

-0.16  -0.04  1.00 

Oak  Grove 

TEMP 

Relhum 

pBARO  SOLRAD  Max  03 

Isoprene 

a-Pinene  3-Pinene  Limonene 

TEMP 

1.00 

RELHUM 

0.04 

1.00 

pBARO 

-0.39 

-0.29 

1.00 

SOLRAD 

0.38 

-0.79 

0.08 

1.00 

Max  Ozone 

0.44 

-0.47 

-0.13 

0.61 

1.00 

Isoprene 

0.67 

0.17 

-0.15 

0.13 

0.12 

1.00 

a-Pinene 

0.45 

0.31 

-0.34 

-0.01 

0.08 

0.63 

1.00 

b-Pinene 

0.57 

0.28 

-0.23 

0.06 

0.21 

0.50 

0.65  1.00 

Limonene 

0.02 

0.22 

-0.03 

-0.16 

-0.21 

0.02 

0.00  0.07  1 

Yorkville 

TEMP 

Relhum 

pBARO  SOLRAD  Max  03 

Isoprene 

a-Pinene  b-Pinene  Limonene 

TEMP 

1.00 

RELHUM 

-0.07 

1.00 

pBARO 

-0.12 

-0.30 

1.00 

SOLRAD 

0.33 

-0.65 

0.26 

1.00 

max  ozone 

0.81 

-0.26 

-0.06 

0.41 

1.00 

Isoprene 

0.76 

0.00 

0.06 

0.18 

0.63 

1.00 

a-Pinene 

0.17 

0.49 

-0.23 

-0.54 

-0.06 

0.17 

1.00 

b-Pinene 

0.60 

0.31 

-0.11 

-0.05 

0.55 

0.59 

0.38  1.00 

Limonene 

-0.19 

0.23 

0.12 

-0.39 

-0.21 

-0.09 

0.40  0.12  1 

Candor 

TEMP 

Relhum 

Max  03 

Solar 

Isoprene 

a-Pinene 

b-Pinene  Limonene 

TEMP 

1.00 

RELHUM 

-0.04 

1.00 

max  ozone 

0.63 

-0.51 

1.00 

SOLRAD 

0.31 

-0.71 

0.62 

1.00 

Isoprene 

0.62 

-0.02 

0.43 

0.25 

1.00 

Benzene 

0.38 

0.25 

0.00 

-0.07 

0.16 

a-Pinene 

0.43 

0.53 

-0.03 

-0.26 

0.34 

1.00 

b-Pinene 

0.31 

0.55 

-0.02 

-0.45 

0.10 

0.70 

1.00 

Limonene 

-0.34 

0.28 

-0.37 

-0.39 

-0.03 

0.36 

0.30  1.00 
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3.3.  Effect  of  Hydrocarbon  Reactivity  with  Hydroxyl  Radical 

Characterizing  the  abundance  of  hydrocarbons  at  a  site  provides  information 

concerning  source  contribution,  however,  it  does  not  take  into  account  individual 
compound  reactivities.  The  contribution  of  the  hydrocarbon  compounds  to  the  production 
of  photochemical  ozone  is  related  to  their  reaction  with  hydroxyl  radicals  and  ozone  in  the 
complex  photooxidation  mechanism.  For  most  hydrocarbons,  the  OH  radical  is  the  most 
important  reaction  pathway.  While  a  species  may  have  a  high  concentration  at  a  given 
site,  if  it  is  not  highly  reactive  with  the  hydroxyl  radical,  then  it  does  not  play  as  important 
a  role  as  a  precursor  to  ozone  formation  as  a  compound  with  a  high  reaction  rate.  To 
account  for  the  reactivity  of  the  hydrocarbon  as  well  as  its  concentration,  we  have  adopted 
the  method  used  by  Chamedies  et  al.,  1992,  and  Lawrimore  et  al.,  1995,  by  calculating  the 
propylene-equivalent  concentration: 


k 

Propy-Equiv (j)  =  Conc(j)  ,  . 


Propy-Equiv(j)  is  a  measure  of  the  concentration  of  species  j  on  an  OH-reactivity  based 
scale,  normalized  to  the  reactivity  of  propylene,  Conc(j)  is  the  concentration  of  species  j 
in  ppb  of  carbon,  ^ohO)  is  the  rate  constant  for  the  reaction  between  species  j  and  OH,  and 
^ohCCsHc)  is  the  rate  constant  for  the  reaction  between  OH  and  propylene  [Chameides  et 
al.,  1992].  The  propylene-equivalent  concentration  is  literally  the  concentration,  in  ppbC, 
required  of  propylene  to  yield  a  carbon  oxidation  rate  equal  to  that  of  the  species].  For 
example,  if  a  species  j  had  a  concentration  of  5  ppbC  and  was  twice  as  reactive  as 
propylene,  it  would  have  a  propylene-equivalent  concentration  of  10  ppbC.  This  method 
is  usefiil  since  it  accounts  for  the  reaction  rate  of  a  species  as  well  as  its  atmospheric 
concentration.  Table  6  lists  the  rate  constant  k  multiplied  by  10*^  for  the  gas-phase 
reactions  of  the  OH  radical  with  hydrocarbons.  Units  of  k  are  in  cm^  molecule'*  s'*.  The 
OH  rate  constants  for  some  of  the  hydrocarbon  compounds  could  not  be  found  in  the 
literature  and  were  therefore  omitted  from  the  propylene-equivalent  hydrocarbon  sums. 
These  compounds  include  3 -methyl- 1-pentene,  4-methy-l-pentene,  cis-3-hexene. 
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Table  6.  Rate  constants  k  (cm^  molecule'^  s’*)  for  the  reactions  of  OH  radicals  with 
hydrocarbons  at  T  =  298°K  (from  Atkinson,  1990,  except  where  noted). 


Compound 

10'^  X* 

Compound 

lO'^xA: 

Ethene 

8.52 

2,4-Dimethylpentane 

5.1 

Acetylene 

0.9 

1, 1, 1-Trichloroethane^ 

0.0119 

Ethane 

0.268 

Benzene 

1.23 

Propene 

26.3 

Cyclohexane 

7.49 

Propane 

1.15 

2,3-Dimethylpentane 

Isobutane* 

2.5 

Trichloroethylene* 

2.36 

Isobutene* 

51 

Methylcyclohexane 

10.4 

1-Butene* 

31.4 

2,3,4-Trimethylpentane 

7 

n-Butane 

2.54 

Toluene 

5.96 

Trans-2-butene 

64 

n-Octane 

8.68 

Cis-2-butene 

56.4 

Perchloroethylene* 

2.16 

3-methyl-l-Butene 

31.8 

Ethylbenzene 

7.1 

Isopentane* 

3.1 

p-Xylene 

14.3 

1-Pentene 

31.4 

m-Xylene 

23.6 

2-methyl- 1  -Butene 

61 

Styrene* 

5.71 

n-Pentane 

3.94 

o-Xylene 

13.7 

Isoprene 

101 

Isopropylbenzene 

6.5 

Trans-2-pentene 

67 

a-Pinene 

53.7 

Cis-2-pentene 

65 

n-Propylbenzene 

6 

2-methyl-2-Butene 

68.9 

l-ethyl-3-Methylbenzene* 

22.4 

3  -methyl- 1  -Pentene 

l-ethyl4-Methylbenzene* 

13.6 

4-methyl- 1  -Pentene 

1,3,5-Trimethylbenzene 

57.5 

Cyclopentane 

5.16 

l-ethyl-2-Methylbenzene* 

13.2 

2-Methylpentane 

5.6 

b-Pinene 

78.9 

3  -Methylpentane 

5.7 

1,2,4-Trimethylbenzene* 

37.23 

n-Hexane 

5.61 

Limonene 

170 

Cis-3 -hexene 

1,3-Diethylbenzene 

Methylcyclopentane^ 

6.596 

n-Butylbenzene 

*  Wameck  1988 
^  Middleton  and  Stockwell,  1990 

^NIST  Chemical  Kinetics  Database,  Version  5.0  [Westley  et  al.,  1993] 


2,3-dimethylpentane,  1,3-diethylbenzene  and  n-butylbenzene.  Because  these  compounds 
were  left  out,  the  propylene-equivalent  concentrations  of  the  sums  may  be  somewhat 
underestimated,  though  not  by  much  since  the  sum  of  these  6  compounds  range  between  ~ 
3  and  4.5  ppbC  at  the  four  sites  during  the  summer  and  between  -1.5  and  2  ppbC  during 
the  winter. 

Figure  7  a-b  shows  selected  hydrocarbons  in  propylene-equivalent  concentrations. 
The  hydrocarbon  sums  were  calculated  by  first  individually  calculating  the  propylene- 
equivalent  concentration  for  each  hydrocarbon,  then  summing  them  up  in  their  appropriate 
categories.  Figure  7a  shows  that,  by  taking  reactivity  into  account,  the  summer 
contribution  of  the  biogenics  at  Candor  is  at  least  65%  of  the  total  sum  of  hydrocarbons 
given  as  propy-equivalent  concentrations.  The  biogenics  include  isoprene,  a-pinene,  P- 
pinene,  and  limonene.  The  biogenic  contribution  was  highest  at  the  Centreville  site, 
contributing  90%  to  the  total  sum.  The  biogenics  at  the  Oak  Grove  and  Yorkville  sites 
represented  78%  and  69%  of  the  total,  respectively.  Isoprene  was  clearly  the  dominant 
compound  during  the  summer,  having  a  propylene-equivalent  concentration  ranging  from 
81.22  ppbC  at  the  Centreville  site,  to  37.63  ppbC  at  the  Yorkville  site.  The  reaction 
between  isoprene  and  the  hydroxyl  radical  is  approximately  3.84  times  faster  than  the 
reaction  between  propylene  and  OH.  This  means  that  at  the  Centreville  site,  for  example, 
a  concentration  of  81  ppbC  of  propylene  would  be  required  to  yield  a  carbon  oxidation 
rate  equivalent  to  21  ppbC  of  isoprene.  Also  from  figure  7a,  the  propy-equivalent 
concentrations  for  the  paraffins,  olefins  and  aromatics  are  much  lower  than  the  biogenic 
compounds.  For  example,  the  propylene-equivalent  alkane  sums  range  from  2  ppbC  to 
5.5  ppbC  among  the  sites.  Table  7  includes  the  ten  most  abundant  species  in  propylene- 
equivalent  concentrations  during  the  summer  and  winter  seasons  for  each  site. 

In  contrast  to  the  summer  propylene-equivalent  concentrations,  the  total 
propylene-equivalent  nonmethane  hydrocarbon  sum  was  four  to  seven  times  lower  for  the 
winter  season,  ranging  from  15  to  21  ppbC,  as  shown  in  figure  7b.  Despite  virtually 
negligible  isoprene  emissions  during  the  winter  season,  the  biogenics  still  dominate  in  the 
winter  due  to  the  emissions  of  terpenes  throughout  the  year.  When  taking  reactivity  into 
account,  limonene  was  the  dominant  terpene  compound  at  all  sites,  ranging  from  -3  to  -  4 
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Propy-Eqivalent  Concentrations:  Summer 


Table  7.  Ten  dominant  hydrocarbons  calculated  in  propylene-equivalent  concentrations,  listed  in  descending  order. 
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ppbC.  This  is  in  sharp  contrast  to  the  regular  concentrations,  in  which  the  alkanes 
dominate.  While  the  paraffins  range  from  28  to  35  ppbC  during  the  winter,  their 
calculated  propylene-equivalent  concentrations  are  approximately  3  ppbC.  It  must  be 
pointed  out  that  the  propylene-equivalent  approach,  or  any  other  OH  reactivity  concept, 
prioritizes  the  individual  compounds  in  terms  of  producing  RO2  radicals.  However,  the 
availability  of  NOx  is  essential  for  the  production  of  photochemical  ozone.  Also,  the 
ozone  present  at  these  rural  sites  will  compete  with  OH  for  the  reaction  with  the  olefin 
compounds,  complicating  the  reaction  mechanism  of  ozone  production.  Although  the  rate 
coefficients  for  reactions  between  olefins  and  ozone  are  much  smaller  than  those  between 
olefins  and  the  OH  radical,  the  reactions  become  competitive  with  OH  when  the 
concentration  of  ozone  builds  up  (Wameck  1988,  189). 
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3.4.  Relationship  between  Trace  Gas  and  Meteorological  Variables 

Meteorology  plays  an  important  role  in  the  Southeast’s  ozone  problem.  Stagnant 
high  pressure  systems  often  develop  over  the  southeast  in  the  summer  and  fall,  which 
allows  a  steady  buildup  of  ozone  precursors.  These  stagnating  high  pressure  systems 
often  arise  from  the  stalling  of  a  continental  high-pressure  system  over  the  Appalachian 
Mountains  [SOS  1993,  12].  This  is  possibly  influenced  by  the  semi-permanent  Bermuda 
High,  which  is  normally  located  over  the  Atlantic  Ocean  [Aneja  and  Yoder  1992,  14]. 

The  summer  of  1993  was  particularly  conducive  to  ozone  formation.  June  and 
July  was  dominated  by  a  persistent  circulation  pattern  which  brought  moisture  from  the 
Ciulf  region  to  the  Midwest,  resulting  in  excessive  precipitation  which  caused  severe 
flooding,  while  preventing  the  eastward  progression  of  weather  systems  which  would  have 
brought  rain  and  cooler  weather  to  relieve  the  Southeast  U.S.  of  drought  conditions  and 
record  breaking  high  temperatures.  The  Southeast  experienced  the  second  warmest 
summer  in  99  years  of  record  keeping*.  For  the  week  of  July  18-24,  above  normal 
temperatures  affected  the  area  from  the  southeastern  Plains  to  the  middle  and  southern 
Atlantic  Coast,  with  temperatures  averaging  3°  to  8°  higher  than  normal.  Abnormally  hot 
weather  continued  through  July  25-3 1^.  Figures  8a-g  show  surface  and  500  mb  height 
contours  for  July  19th  through  the  25th.  Throughout  this  period,  an  upper  level  high 
characterized  by  light  500  mb  winds  was  centered  over  the  Southeast.  During  the  period 
of  the  ozone  episode  (July  20-24),  a  stationary  blocking  pattern  characterized  by  a  high 
amplitude  ridge  or  an  omega  block  was  centered  over  the  central  states,  then  over  the 
Mississippi  valley  region.  An  omega  block  is  a  high  amplitude  ridge  shaped  as  the  capital 


‘  U.S.  Dept,  of  Commerce.  NOAA.  National  Weather  Service.  Climate  Analysis  Center.  United  States  seasonal 
climate  summary,  summer  (June-August)  1993.  Weekly  Climate  Bulletin.  Ed.  Richard  Tinker.  No.  93/39. 
Washington  DC:  Govt.  Printing  Office,  September  29,  1993. 

^  U.S.  Dept,  of  Conunerce.  NOAA.  National  Weather  Service.  Climate  Analysis  Center.  United  States  weekly 
climate  highlights,  for  the  week  of  July  18-24,  1993.  Weekly  Climate  Bulletin.  Ed.  Richard  Tinker.  No.  93/30. 
Washington  DC:  Govt.  Printing  Office,  July  28,  1993. 

^  U.S.  Dept,  of  Commerce.  NOAA.  National  Weather  Service.  Climate  Analysis  Center.  United  States  weekly 
climate  highlights,  for  the  week  of  July  25-31,  1993.  Weekly  Climate  Bulletin.  Ed.  Richard  Tinker.  No.  93/31. 
Washington  DC:  Govt.  Printing  Office,  August  4,  1993. 
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Figures  8  d-e.  500  mb  height  contours  at  7am  EST  on  (d)  Thursday  July  22,  1993  (top 
figure)  and  (e)  Friday  July  23,  1993  (bottom  figure). 


500  mb  height  contours  at  7am  EST,  Sunday  July  25,  1993 


Greek  letter  n,  and  is  sandwiched  in  between  two  lows  [Bluestein  1993,  79],  In  a 
blocking  pattern,  those  living  near  upper-level  cyclones  tend  to  experience  a  persistent 
combination  of  precipitation  and  relatively  cool  temperatures,  while  those  near  upper-level 
anticyclones  tend  to  experience  drought  conditions,  which  was  the  situation  for  this 
period. 

Figure  9  is  a  plot  of  hourly  averaged  ozone  during  an  episode  which  occurred  from 
July  20th  through  July  25th,  1993.  Ozone  values  for  July  20th  to  22nd  at  the  Candor  site 
are  not  included  because  of  system  malfunction.  On  July  19,  a  warm  front  passed  through 
the  southeastern  Atlantic  coast  states.  The  7  am  EST  surface  map  shows  the  warm  front 
situated  across  Virginia,  North  and  South  Carolina.  Ozone  was  low  at  all  of  the  sites  on 
this  day  (Figure  9).  Yorkville  received  3.8  cm  of  precipitation  the  previous  evening  on  the 
18th,  and  4  cm  the  next  morning  on  the  19th  as  a  result  of  the  frontal  passage.  Candor 
received  precipitation  later  in  the  day  from  10  am  to  4  pm.  Oak  Grove  also  received 
precipitation  late  in  the  afternoon,  from  4  to  5  pm.  Centreville  did  not  receive  any  rain 
during  the  episode.  By  comparing  the  diurnal  profiles  of  Centreville  and  Yorkville,  it 
appears  that  ozone  at  Centreville  preceded  that  of  Yorkville  by  approximately  a  day,  its 
maximum  increasing  a  day  earlier  and  declining  two  days  sooner  than  Yorkville.  On  the 
20th,  the  high  pressure  center  over  the  southeastern  states,  influenced  by  the  omega  block 
over  northern  Canada  and  upper  level  ridge  over  the  mid-west,  is  the  only  distinct  weather 
feature  which  would  explain  the  increase  in  ozone  concentrations  at  the  Centreville  and 
Yorkville  sites.  On  the  21st,  the  omega  block  had  moved  from  northern  Canada  down 
towards  the  central  plains  of  the  U.S.,  while  the  high  pressure  and  high  temperatures 
continued.  Daily  maximum  ozone  exceeded  100  at  both  Centreville  (109  ppb)  and 
Yorkville  (108  ppb)  while  increasing  to  54  ppb  at  Oak  Grove.  On  the  22nd,  the  omega 
block  became  quite  distinctive,  having  moved  towards  the  east  and  centering  over  the 
Mississippi  Valley.  High  pressure  continued  over  the  southeast  with  light  500  mb  winds 
and  high  temperatures.  Daily  maximum  ozone  exceeded  100  ppb  at  both  Centreville  and 
Yorkville  for  a  second  day,  reaching  104.5  ppb  at  Centreville  and  peaking  at  1 15  ppb  at 
Yorkville.  Ozone  reached  65  ppb  at  Oak  Grove.  On  the  23rd  and  24th,  the  upper  level 
omega  blocking  pattern  remained  centered  over  the  Mississippi  Valley  region.  High 
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Figure  9.  Diurnal  profiles  of  ozone  during  an  episode  which  occurred  July  20-25  1993. 


pressure  and  high  temperatures  continued  in  the  Southeast.  Daily  maximum  ozone  on  the 
23rd  decreased  to  84.5  ppb  at  Centreville  and  to  58.75  at  Oak  Grove,  while  peaking  at 
1 10.5  ppb  at  Yorkville.  Ozone  decreased  at  Centreville  and  Oak  Grove  on  July  24  while 
remaining  high  at  Yorkville  (1 13  ppb  maximum).  On  July  25,  a  trof  line  shown  on  the 
NWS  surface  analysis  shows  it  passing  through  the  southeastern  Atlantic  coast  states. 

Daily  maximum  ozone  reached  only  65  ppb  at  Yorkville,  which  can  be  explained  by 
precipitation  which  occurred  throughout  the  afternoon  and  evening  on  the  previous  day 
(July  24)  and  again  through  the  afternoon  and  evening  on  the  25th.  Ozone  remained  low 
at  Centreville  (49  ppb). 

Figure  10  shows  the  trace  gases  (O3,  CO,  SO2,  NOy)  measured  at  the  sites  during 
the  hydrocarbon  sampling  period  (every  6  days,  1200  to  1300).  The  only  trace  gas 
measured  at  Candor  was  ozone.  From  figure  10  we  see  that  the  Yorkville  site  has  the 
highest  summer  average  of  ozone  of  all  the  sites  (93  ±  22  ppb).  It  also  has  the  highest 
CO,  SO2,  NO,  and  NOy  among  the  three  sites  for  which  measurements  were  taken.  This 
indicates  that  Yorkville,  while  perhaps  located  in  a  rural  site,  is  subject  to  the  influence  of 
anthropogenic  emissions.  Yorkville  is  located  approximately  45  miles  (72  km)  west  of 
Atlanta,  Georgia,  and  pollutant  transport  is  most  likely  the  reason  for  these  high  values. 
The  average  ozone  for  Centreville  and  Oak  Grove  during  the  summer  remained  essentially 
the  same  as  during  the  spring  and  autumn. 

Table  8  shows  the  mean  and  standard  deviation  of  the  meteorological  and  trace  gas 
data  for  the  period  from  June  1  to  August  3 1,  1993,  using  data  averaged  from  10  am  to  4 
pm  for  each  day,  using  approximately  92  observations.  As  can  be  seen  from  Table  8, 
ozone  was  the  highest  at  the  Yorkville  and  Candor  sites.  Relative  humidity,  CO,  SO2, 

NO,  and  NOy  were  also  highest  at  Yorkville.  A  t-test  was  performed  for  these  parameters 
between  Yorkville  and  Centreville,  and  Yorkville  and  Oak  Grove.  The  t-test  showed  that 
daily  average  ozone  at  Yorkville  was  significantly  higher  than  at  Centreville  or  Oak  Grove 
at  the  1%  level  (p  =  0.0001).  SO2,  NO,  and  CO  also  showed  up  as  significantly  higher  at 
the  Yorkville  site  (p=0.0001).  The  difference  in  NOy  between  Yorkville  and  Oak  Grove 
was  significant  at  the  1%  level,  and  significant  at  the  5%  level  between  Yorkville  and 
Centreville.  The  higher  concentrations  of  these  compounds  shows  that  Yorkville,  despite 
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Figxire  10.  Seasonal  averages  of  trace  gases  measured  during  the  hydrocarbon  sampling  period.  Ozone  represents 
the  seasonal  average  of  daily  maximum  ozone. 


Table  8.  Descriptive  statistics  of  meteorological  and  trace  gas  data  for  all  four  sites.  The 
data  used  represents  daily  averages  from  1000  to  1600,  for  the  period  June  1  through 
August  31,  1993  (n  =  92).  The  daily  max  O3  represents  the  maximum  ozone 


concentration  occurring  between  1000  and  1600  each  day. 


Summer  1993 

Centreville 

Mean  StDev 

Oak  Grove 

Mean  StDev 

Yorkville 

Mean 

StDev 

Candor 

Mean 

StDev 

Temp  (“C) 

30.58 

2.51 

29.72 

2.38 

28.01 

2.62 

29.30 

5.62 

Relative  Humidify 

62.55 

11.93 

64.30 

11.98 

67.51 

11.68 

44.59 

10.31 

pBARO  (mmHg) 

746.14 

2.43 

754.65 

1.46 

728.95 

1.74 

Solar  Rad  (w/m^) 

585.96 

176.45 

577.07 

197.43 

599.36 

262.13 

592.12 

256.55 

daily  avg  O3  (ppb) 

47.57 

14.66 

43.63 

11.92 

63.47 

16.58 

65.00 

11.54 

daily  max  O3  (ppb) 

56.19 

17.88 

51.22 

13.54 

76.13 

22.67 

CO  (ppb) 

164.70 

45.84 

184.31 

50.24 

297.68 

58.25 

SO2  (ppb) 

2.11 

2.58 

1.95 

1.84 

5.39 

5.94 

NO  (ppb) 

0.13 

0.13 

0.13 

0.21 

0.30 

0.25 

NOy  (ppb) 

4.47 

2.48 

2.58 

0.98 

5.42 

2.41 

Table  9.  Correlations  between  meteorological  variables  and  trace  gases.  The  data  used 
are  daily  averages  for  the  period  June  1  through  August  3 1,  1993.  Avg  O3  is  the  ozone 
concentration  averaged  from  1000  to  1600  each  day,  and  max  O3  is  the  daily  maximum 


ozone  concentration  occurring  between  1000  and  1600. 


York  93 

Temp 

RH 

pBaro 

Solar 

avgOs 

max  O3 

CO 

SO2 

NO 

NOy 

TEMP 

1.00 

RELHUM 

-0.75 

1.00 

pBARO 

-0.19 

0.05 

1.00 

SOLRAD 

0.31 

-0.51 

-0.09 

1.00 

avg03 

0.51 

-0.44 

0.03 

0.00 

1.00 

max03 

0.48 

-0.41 

0.07 

-0.07 

0.95 

1.00 

CO 

-0.04 

0.18 

-0.03 

0.21 

0.38 

0.35 

1.00 

SO2 

0.17 

-0.20 

0.12 

0.02 

0.36 

0.40 

0.07 

1.00 

NO 

-0.14 

0.05 

0.23 

-0.01 

-0.10 

-0.02 

0.13 

0.36 

1.00 

NOy 

0.26 

-0.22 

0.16 

-0.14 

0.75 

0.78 

0.47 

0.53 

0.31 

1.00 
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Figure  11.  Ozone  versus  NOy-NO.  Data  used  were  daily  averages  from  1000  to  1600,  June  1  through  August  31  1993  (92  observations) 


its  rural  location,  is  heavily  impacted  by  anthropogenic  pollutant  sources,  probably  as  a 
result  of  its  relative  proximity  to  Atlanta,  Georgia,  and  to  the  power  generating  station  to 
the  east.  T-tests  on  the  meteorological  variables  show  that  temperature  and  pressure  were 
significantly  lower  at  Yorkville  than  at  the  other  two  sites  (p=0.0001).  Solar  radiation 
was  not  found  to  be  statistically  different. 

Daily  average  ozone  was  plotted  against  the  difference  NOy-NO  in  Figure  11. 
Ideally,  one  would  plot  the  difference  NOy-NO*,  where  NO*  =  NO  +  NO2,  because  this  is 
a  direct  measure  of  the  products  of  the  NO*  oxidation  and  minimizes  the  variability  due  to 
differences  in  photochemical  aging  of  the  sampled  air  mass  [Trainer  et  al.,  1993]. 
Unfortunately,  NO2  was  not  measured  and  NOy-NO*  could  not  be  plotted,  and  therefore 
the  age  of  the  air  mass  was  not  taken  into  account.  A  correlation  can  be  seen  between 
ozone  and  NOy-NO  at  each  site.  Observed  ozone  at  Yorkville  is  higher  than  that  observed 
at  Centreville  for  a  given  NOy-NO  value.  This  may  be  explained  by  the  fact  that  Yorkville 
is  45  miles  west  of  a  large  anthropogenic  area  source  (Atlanta,  Georgia),  and  30  miles 
west  of  a  large  power  generating  station,  and  is  affected  by  air  containing  relatively 
unaged  NOy,  which,  in  the  presence  of  biogenic  hydrocarbons,  can  result  in  high  ozone 
concentrations.  It  can  also  be  seen  that  Oak  Grove  is  a  particularly  clean  site,  with  daily 
average  ozone  values  not  exceeding  ~  70  ppb  and  NOy-NO  not  exceeding  ~  5  ppb,  so  the 
regression  line  for  Oak  Grove  is  much  steeper  and  has  a  lower  intercept.  Earlier  in  this 
work  it  was  shown  that  during  the  summer,  the  Yorkville  site  had  lower  concentrations  of 
isoprene  (Figure  5  a)  and  non-methane  hydrocarbons  calculated  in  propylene-equivalent 
concentrations  (Figure  7a)  than  at  Centreville  or  Oak  Grove,  yet  Yorkville  had  the  highest 
concentrations  of  ozone.  Ozone  production  in  rural  areas  is  typically  limited  by  the 
availability  of  NO*  rather  than  hydrocarbons,  since  isoprene  and  other  biogenic  VOCs 
provide  a  ubiquitous  source  of  hydrocarbon  precursors  for  ozone  production.  The  high 
ozone  concentrations  at  Yorkville  highlights  this  NO*  limited  characteristic,  showing  that 
greater  concentrations  of  reactive  nitrogen  (NOy)  at  Yorkville  play  a  more  significant  role 
in  ozone  formation  than  greater  concentrations  of  highly  reactive  biogenic  hydrocarbons 
(Centreville).  Correlations  between  meteorological  variables  and  trace  gases  at  Yorkville 
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(Table  9)  show  that  there  is  a  strong  correlation  between  NOy  and  daily  averaged  ozone 
(0.75),  and  daily  maximum  ozone  (0.78). 

A  multiple  linear  regression  analysis  was  performed  for  the  summer  1993  ozone 
values  at  Yorkville.  The  variables  which  were  examined  for  the  regression  analysis 
included  hourly  averaged  values  of  temperature,  relative  humidity,  barometric  pressure, 
solar  radiation,  CO,  SO2,  and  NOy.  The  multiple  regression  showed  that  temperature, 
relative  humidity,  and  NOy  together  accounted  for  most  of  the  variation  in  daily  average 
and  daily  maximum  ozone  at  a  1%  significance  level.  For  daily  average  ozone,  the 
regression  equation  becomes 

Ozone  =  1.698(Temperature)  -  0. 13(Relative  humidity)  +  4.55(NOy) 

A  plot  of  the  observed  ozone  values  and  predicted  values  based  on  the  above  equation  is 
shown  in  Figure  12.  As  you  can  see,  the  estimated  ozone  values  follow  the  actual  values 
rather  well,  the  only  exception  to  this  being  the  erroneously  high  prediction  on  August 
1 1th,  which  is  due  to  a  daily  averaged  NOy  spike  of  14.26  ppb.  The  R-squared  value  of 
this  regression  was  reported  in  the  output  of  the  SAS  statistical  program  to  be  0.978, 
though  it  must  be  noted  that  this  R-squared  value  is  misleading.  The  intercept  was 
removed  from  the  equation  because  it  was  not  statistically  significant  at  the  5%  level, 
however,  the  absense  of  the  intercept  term  causes  the  total  sum  of  squares  to  become 
extremely  large  (i.e. :  SS(total)  =  2:(yi  -  y  )^  becomes  SS(total)  =  Z(yi)^  )  and  the  R- 
squared  value  approximates  to  1  (  R^  =  1  -  [SS(error)/SS(total)]  ).  Therefore,  for  the 
purpose  of  obtaining  an  accurate  R-squared  value,  the  multiple  regression  was  re¬ 
evaluated  with  the  intercept  term  included  in  the  model.  This  yielded  an  R-squared  value 
of  0.66,  though  the  intercept  was  still  not  statistically  significant,  and  its  inclusion  in  the 
model  resulted  in  the  relative  humidity  parameter  being  not  significant  at  the  5%  level. 
There  was  virtually  no  change  in  the  predicted  ozone  output  between  the  two  regression 
equations,  so  R  »  0.66  is  used  here  as  an  approximate  R-squared  value  for  the  regression 
equation  reported  above. 
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e  regression  of  ozone  at  YorkvilJe, 


Because  Yorkville  has  greater  concentrations  of  anthropogenic  compounds  than 
Oak  Grove,  a  separate  multiple  linear  regression  was  done  for  the  Oak  Grove  data  to 
compare  parameter  estimates.  Temperature,  relative  humidity  and  NOy  were  found  highly 
significant  in  the  linear  regression  of  daily  average  ozone  at  Oak  Grove  as  well.  For  both 
the  Yorkville  and  Oak  Grove  regressions,  temperature  and  NOy  were  highly  significant  (p 
=  0.0001)  while  relative  humidity  was  significant  at  the  5%  level  (p  =  0.033  for  Oak 
Grove  and  0.037  for  Yorkville).  The  regression  equation  for  daily  average  ozone  at  Oak 
Grove  was 

Ozone  =  1.181(Temperature)  -  0.134(Relative  humidity)  +  6.497(NOy) 
which  is  very  similar  to  the  Yorkville  regression  equation. 
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Summary 

An  analysis  of  hydrocarbons  sampled  from  1992  through  1993  at  four  rural  sites  in 
the  Southeast  shows  a  seasonal  variation  of  light  molecular  weight  (C2-C4)  alkanes,  ethene 
and  acetylene,  with  a  maximum  during  the  winter  and  minimum  during  the  summer.  The 
biogenic  hydrocarbons  (isoprene  and  the  terpenes)  also  display  a  seasonal  variation,  with  a 
summer  maximum  and  winter  minimum.  Isoprene  was  virtually  non-existent  during  the 
winter  at  all  sites,  and  averaged  from  9.8  ppbC  (Yorlcville,  GA)  to  21.15  ppbC 
(Centreville,  AL)  during  the  summer.  The  terpene  concentration  was  greatest  in  the 
summer  with  averages  ranging  3.19  ppbC  (Centreville,  AL)  to  6.38  ppbC  (Oak  Grove, 
MS),  but  was  also  emitted  during  the  winter  months,  with  a  range  of  1.25  to  1.9  ppbC  for 
all  sites.  When  considering  the  reactivity  of  hydrocarbons  with  the  OH  radical,  the 
biogenics  dominate  the  total  non-methane  hydroarbon  sum,  representing  between  65%  to 
90%  of  TNMHCs  during  the  summer  season,  while  the  impact  of  the  other  hydrocarbons 
are  less  important.  The  propy-equivalent  TNMHC  sums  during  the  summer  at  the  four 
sites  range  between  70  to  100  ppbC  with  isoprene  being  the  dominant  hydrocarbon;  this 
propy-equivalent  range  drops  during  the  winter  season  with  a  range  of  15  to  20  ppbC, 
when  isoprene  emissions  are  negligible. 

Seasonal  averages  of  the  trace  gases  show  that  Yorkville  was  the  most  affected  by 
anthropogenic  emissions,  while  Oak  Grove  was  the  cleanest  of  the  sites.  Despite  the  fact 
that  Yorkville  had  the  lowest  concentration  of  summer  propy-equivalent  total  NMHCs,  it 
had  the  highest  values  of  ozone,  SO2,  NOy  and  CO  than  the  other  two  rural  sites  for  which 
measurements  were  taken.  A  plot  of  ozone  versus  NOy-NO  shows  that  with  a  given 
concentration  of  NOy-NO,  the  Yorkville  site  had  higher  ozone  concentrations  than  either 
the  Centreville  or  Oak  Grove  sites.  These  observations  highlight  the  NOx  limited 
characteristic  of  this  region,  namely,  that  because  hydrocarbons  are  ubiquitous  in  this 
region  due  to  natural  hydrocarbon  emissions,  the  ozone  producing  potential  is  limited  to 
the  availability  of  NOx,:  in  this  region,  greater  levels  of  NO*  play  a  more  significant  role 
inozone  formation  than  greater  concentrations  of  reactive  biogenic  hydrocarbons. 
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SECTION  n 


MULTIPLE  REGRESSION  OF  OZONE  USING  TRACE  GAS  AND 
METEOROLOGICAL  VARIABLES  WITH  AUTOREGRESSIVE 
MOVING  AVERAGE  (ARMA)  TIME  SERIES  ERRORS 


1.  Introduction 

Ozone  (O3)  is  an  oxidant  gas  produced  naturally  in  the  atmosphere.  In  the 
troposphere,  it  is  produced  when  nitrogen  dioxide  (NO2)  is  photodisassociated  and 
combines  with  molecular  oxygen  (O2)  to  produce  ozone.  Since  the  passage  of  the  1970 
Clean  Air  Act  amendments,  regulatory  efforts  to  comply  with  the  National  Ambient  Air 
Quality  standard  for  ozone  have  been  inadequate  [NRC  1991,4].  Ozone  exceedences 
continue  to  be  a  major  problem,  especially  in  the  southeast  region  of  the  United  States 
which  is  characterized  by  high  frequency  of  stagnating  high  pressure  systems  and 
significant  emissions  of  naturally  produced  hydrocarbons  (SOS  1995,  12). 

Many  studies  have  attempted  to  model  the  formation  of  ozone  using  highly 
complex  computer  simulations  of  chemical  reactions  of  ozone  precursors  [Chameides  et 
al.,  1988;  Trainer  et  al.,  1987].  Highly  involved  models  incorporating  atmospheric 
chemistry  as  well  as  meteorological  variables  are  necessary  for  accuracy.  However,  it  is 
interesting  to  determine  if  a  statistical  model  using  observed  ozone  concentrations  and 
meteorological  input  parameters  at  a  given  site  can  simulate  ozone  concentrations  for 
another  year  with  moderate  accuracy. 

In  this  work,  a  multiple  regression  using  meteorological  and  trace  gas  data  as  input 
parameters  and  ARMA  (autoregressive  moving  average)  time  series  errors  was  applied  to 
model  daily  average  ozone.  The  data  used  was  collected  from  June  1  through  August  3 1 
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1992  at  Yorkville,  a  rural  site  located  in  Pauldin  County,  Georgia  (33°55’4r’  N, 
85°02’46”  W).  This  statistical  model  was  then  applied  to  the  same  site  for  the  1993 
summer,  and  the  predicted  and  observed  ozone  values  were  compared  using  statistical 
tests. 


2.  Time  Series  Review:  Autoregressive  moving  average  (ARMA)  processes 

Autocorrelation  is  the  correlation  of  a  variable  with  itself  between  the  variable 
sampled  at  time  t  and  later  at  time  t  +  lag  (the  correlation  between  Yt  and  Yt+k).  When 
autocorrelation  is  zero,  Yt  is  a  random  process,  while  perfect  correlation  is  denoted  by  a 
+1  or  -1.  The  correlation  between  Yt  and  Yt+k  is  calculated  by  the  autocorrelation 
function: 


Cov(Y,Y.^,)  _  Yk 

VVar(Y,)VVar(Y,,J  Y„ 

where  Yo  =  Var(Yt)  =  Var(Yt+k) 

Yk  is  called  the  autocovariance  function,  and  pk  is  called  the  autocorrelation  function 
(ACF).  They  represent  the  covariance  and  correlation  between  Yt  and  Yt+k  from  the  same 
process  separated  only  by  a  time  lag  of  k. 

A  time  series  process  can  be  expressed  in  two  ways:  as  an  autoregressive 
representation  (AR),  and  as  a  moving  average  (MA)  representation.  An  autoregressive 
(AR)  process  is  when  a  series  Yt  can  be  predicted  using  past  values  (Yu,  Yt.2,...)  plus  a 
random  shock  (et),  such  as  a  weather  forecast  for  tomorrow  can  be  predicted  based  on  its 
pattern  today  and  yesterday,  with  a  given  amount  of  uncertainty.  A  first  order 
autoregressive  process  includes  one  past  value  of  Yt,  and  the  following  equation  is  an 
example  of  an  AR(1)  process: 


Yt  =  ttiYt-i  +  et 


(2) 
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where  the  ai  coefficient  is  the  weight  of  the  past  value,  and  et  is  the  random  error. 

A  moving  average  (MA)  representation  is  when  the  process  Yt  is  represented  as  a 
linear  combination  of  a  sequence  of  uncorrelated  random  variables,  taking  a  weighted 
average  of  the  errors  where  past  errors  are  less  important.  A  first  order  MA  process 
includes  one  lagged  error  term: 

Yt=  et-pie,.i  (3) 

The  ARMA  (p,q)  model  is  a  mix  of  the  AR  and  MA  processes  together  in  the  same 
equation,  where  p  is  the  order  of  the  AR  process  and  q  is  the  order  of  the  MA  process. 
The  ARMA  (1,1)  model  is  written  as 

Yt-aiYn  =  et-Piet.i  (4) 

The  equation  for  a  regression  with  time  series  errors  is; 

Yt  =  cOo  +  ©iXit  +  0)2X21 +  ...+ CDkXkt  +  Zt  (5) 

where  Xu,  X2t,  etc.,  are  input  variables  and  Oi,  (02...  are  the  unknown  coefficients.  Zt  is  an 
ARIMA  time  series.  This  is  a  typical  regression  model,  except  that  it  allows  for 
autocorrelation  in  the  error  term  Z.  In  this  work,  an  ARMA(1,1)  equation  was  used  for 
the  error  term  Z. 


3.  Results  and  Discussion 

The  descriptive  statistics  and  correlation  coefficients  of  the  ozone,  trace  gas  and 
meteorological  data  from  the  Yorkville  site  during  summer  1992  are  displayed  in  tables  1 
and  2,  respectively.  The  data  was  averaged  from  10  a.m.  to  4  p.m.  local  standard  time  for 
each  day,  except  for  max  O3,  which  represents  the  daily  maximum  ozone  concentration. 
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Yorkville  is  affected  by  two  potential  emission  sources;  the  city  of  Atlanta,  located 
approximately  45  miles  southeast  of  the  site,  and  a  large  power  generating  station  located 
approximately  30  miles  east  of  the  site  [Kirk  1996], 


Table  1.  Descriptive  statistics  of  the  Yorkville  1992  data 


Yorkville  1992 

Mean 

Std  Enor 

Median 

StDev 

Count 

Temperature  C° 

24.80 

0.36 

24.86 

3.41 

88 

Relative  Humidity 

74.12 

1.25 

73.94 

12.05 

93 

pBARO 

728.86 

0.29 

729.20 

2.78 

92 

Solar  Radiation 
(watts  m'^) 

539.70 

22.11 

563.26 

182.33 

68 

Ozone  (ppb) 

50.25 

1.55 

49.70 

14.80 

91 

Max  Ozone  (ppb) 

60.07 

1.80 

59.00 

17.07 

90 

CO  (ppb) 

204.42 

5.01 

204.56 

47.77 

91 

S02  (ppb) 

4.98 

0.66 

2.66 

6.29 

92 

NO  (ppb) 

0.85 

0.14 

0.44 

1.31 

93 

NOy  (ppb) 

8.94 

0.60 

7.88 

5.82 

93 

Table  2.  Correlation  coefficients  between  1992  variables.  The  values  in  boldface  indicate 


correlations  greater  than  0.50  or  less  than  -0.60. 


York  92 

TEMP 

RELHUM 

pBARO 

SOLRAD 

OZONE 

maxOs 

CO 

S02 

NO  NOy 

TEMP 

1.00 

RELHUM 

-0.52 

1.00 

pBARO 

0.35 

0.00 

1.00 

SOLRAD 

0.55 

-0.88 

-0.11 

1.00 

OZONE 

0.18 

-0.67 

-0.14 

0.54 

1.00 

max  O3 

0.15 

-0.60 

-0.10 

0.47 

0.96 

1.00 

CO 

-0.25 

0.17 

-0.16 

-0.33 

0.19 

0.22 

1.00 

S02 

-0.16 

0.00 

0.13 

-0.09 

0.01 

0.11 

-0.06 

1.00 

NO 

-0.25 

0.05 

0.12 

-0.17 

-0.05 

0.05 

0.20 

0.67 

1.00 

NOy 

-0.09 

-0.17 

0.05 

-0.03 

0.45 

0.52 

0.44 

0.44 

0.77  1.00 

NOy  is  the  sum  of  the  reactive  odd  nitrogen  species  (NOy  =  NO  +  NO2  +  organic  nitrates 
+  inorganic  nitrates).  There  appears  to  be  a  strong  negative  correlation  between  ozone 
and  relative  humidity  (-0.67),  a  fair  correlation  between  ozone  and  solar  radiation  (0.54), 
and  with  reactive  nitrogen  NOy  (0.45).  The  relationship  between  ozone  and  meteorology 
is  consistent  with  the  literature.  For  example,  studies  have  been  done  to  determine  long 
term  ozone  trends  by  removing  the  effect  of  meteorology  [Chock  et  al.,  1982;  Kumar  and 
Chock  1984;  Korsog  and  Wolff,  1991].  Ozone  is  produced  from  photochemical  reactions 
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between  volatile  organic  compounds  and  oxides  of  nitrogen  (NO*  =  NO  +  NO2),  and 
therefore  is  influenced  by  the  level  of  NO*  and  solar  radiation.  The  negative  correlation 
between  humidity  and  solar  radiation  is  high  (-0.88). 

The  variables  which  were  examined  for  the  regression  include  temperature,  relative 
humidity,  barometric  pressure,  solar  radiation,  carbon  monoxide  (CO),  sulfur  dioxide 
(SO2),  and  reactive  nitrogen  (NOy).  Nitric  oxide  (NO)  was  not  used  since  NOy  includes 
NO  (NOy  =  NO  +  NO2  +  organic  nitrates  +  inorganic  nitrates).  A  regression  model  with 
ARMA  (1,1)  errors  was  fitted  to  daily  average  ozone  for  the  period  June  1  -  August  3 1 
1992  (refer  to  Appendix  A  for  SAS  programming  and  computer  output).  The  model  is 

Yt  =  -  0.758RHt  +  0.137Pt+  0.577SO2t  +  1.053NOyt  +  Zt  (6) 

where  Zt  =  0.875Zt.i  +  et  -  0.662et.i 

Yt  is  daily  average  ozone,  RH  is  relative  humidity,  P  is  barometric  pressure,  S02  and  NOy 
are  sulfur  dioxide  (SO2),  and  reactive  nitrogen  (NOy),  respectively.  Temperature,  solar 
radiation,  carbon  monoxide  and  an  intercept  term  was  not  found  to  be  statistically 
significant  in  explaining  the  variation  in  ozone  when  taking  into  account  the  other 
variables.  Figure  1  shows  the  observed  ozone  values  plotted  with  the  predicted  ozone 
using  equation  6  above.  The  predicted  ozone  values  follow  the  actual  values  very  well. 

The  observed  ozone  concentrations  was  well  within  the  predicted  95%  confidence 
interval,  and  the  small  values  of  the  chi-squared  statistic  and  high  probability  values 
(>0.05)  in  the  autocorrelation  check  of  residuals  and  shows  that  the  residuals  are  white 
noise,  indicating  that  the  regression  model  fits  well  (refer  to  Appendix  A,  section  II). 

Equation  6  was  next  applied  to  the  1993  Yorkville  data  by  using  Proc  ARIMA  to 
input  the  values  given  fi-om  the  1992  regression  (see  Section  I,  Appendix  B  for 
programming).  The  option  “noest”  was  used  to  prevent  SAS  from  estimating  a  new 
model.  As  can  be  seen  from  Figure  2,  the  actual  daily  average  ozone  was  underpredicted 
most  of  the  time.  Figures  3a  and  3b  show  two  plots  of  residuals.  Figure  3a  is  a  plot  of 
residuals  with  time.  The  residuals  tended  towards  the  positive  side  of  the  zero  axis. 
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having  values  as  high  as  43  ppb.  A  plot  of  residuals  versus  predicted  ozone  in  figure  3  b 
shows  a  linear  relationship.  Ozone  was  underpredicted  for  actual  values  greater  than  80 
ppb.  An  examination  of  the  SAS  output  shows  a  significant  chi-squared  statistic  and  low 
probabilities  (p<0.05)  in  the  autocorrelation  check  of  residuals,  indicating  that  the  model 
does  not  fit  well  and  that  the  residuals  are  not  white  noise  (Appendix  B,  Section  II).  The 
mean  of  the  residuals  (actual  ozone  -  forecast  ozone)  is  4.91,  and  a  t-test  on  the  residual 
mean  shows  a  significant  t-statistic  of  3.68  and  a  probability  of  0.0004,  indicating  that  the 
residuals  are  not  white  noise  (Appendix  B,  Section  HI).  The  model  found  using  the  1992 
data  can  not  be  used  to  predict  ozone  using  the  1993  data.  A  potential  problem  with  using 
regression  analysis  is  that  high  ozone  concentrations  are  consistently  underpredicted  by 
the  statistical  model.  The  least-squares  fitting  procedure  used  in  linear  regression  is 
designed  to  limit  the  overall  mean  square  error,  and  because  high  ozone  values  (greater 
than  100  ppb)  are  extreme  values  and  occur  rarely,  they  might  not  be  important  in 
determining  the  regression  coefficients  (NRC  1991,  61).  Table  3  and  4  lists  the 
descriptive  statistics  and  correlation  coefficients,  respectively,  of  the  1993  Yorkville  data. 


Table  3.  Descriptive  statistics  of  the  Yorkville  1993  data 


Yorkville  1993 

Mean 

Std  Error  Median 

Std  Dev 

Count 

Temperature 

28.01 

0.28 

28.63 

2.62 

90.00 

Relative  humidity  % 

67.51 

1.23 

65.04 

11.68 

90.00 

Pressiue 

728.95 

0.18 

728.81 

1.74 

90.00 

Solar  Radiation 

599.36 

28.60 

689.48 

262.13 

84.00 

Ozone 

63.47 

1.73 

63.62 

16.58 

92.00 

maxo3 

76.13 

2.36 

73.50 

22.67 

92.00 

CO 

297.68 

6.07 

295.40 

58.25 

92.00 

S02 

5.39 

0.62 

3.51 

5.94 

92.00 

NO 

0.30 

0.03 

0.24 

0.25 

92.00 

NOy 

5.42 

0.25 

5.10 

2.41 

92.00 
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Table  4.  Correlation  coefficients  between  1993  variables. 


York  93 

TEMP 

RelHum 

pBARO 

Solar 

Ozone 

max  03 

CO 

S02 

NO 

NOy 

TEMP 

1.00 

RELHUM 

-0.75 

1.00 

pBARO 

-0.19 

0.05 

1.00 

SOLRAD 

0.31 

-0.51 

-0.09 

1.00 

OZONE 

0.51 

-0.44 

0.03 

0.00 

1.00 

maxoS 

0.48 

-0.41 

0.07 

-0.07 

0.95 

1.00 

CO 

-0.04 

0.18 

-0.03 

0.21 

0.38 

0.35 

1.00 

S02 

0.17 

-0.20 

0.12 

0.02 

0.36 

0.40 

0.07 

1.00 

NO 

-0.14 

0.05 

0.23 

-0.01 

-0.10 

-0.02 

0.13 

0.36 

1.00 

0.26 

-0.22 

0.16 

-0.14 

0.75 

0.78 

0.47 

0.53 

0.31 

1.00 

Temperature,  daily  average  ozone,  and  carbon  monoxide  were  found  to  be  significantly 
higher  in  1993  than  in  1992  (p=0.0001).  However,  reactive  nitrogen  (NOy)  was  found  to 
be  significantly  lower  in  1993  than  in  1992  (p=0.0001)  despite  the  higher  ozone  and 
carbon  monoxide  levels  in  1993.  Relative  humidity  was  also  found  to  be  significantly 
lower  in  1993  (p=0.001).  Barometric  pressure,  sulfur  dioxide  (SO2)  and  solar  radiation 
was  not  found  to  be  significantly  different  between  the  two  years.  By  comparing  1992 
correlations  with  1993  correlations,  we  see  that  the  correlation  between  ozone  and  the 
anthropogenic  pollutants  (CO,  SO2,  NOy)  is  higher  in  1993  than  in  1992.  The  correlation 
between  ozone  and  temperature  was  also  higher  in  1993  (0.51  in  1993  vs  0.18  in  1992). 

Variations  in  ozone  levels  are  highly  influenced  by  synoptic  meteorology 
[Vukovich  et  al.,  1977;  Vukovich,  1994;  Niccum  et  al.,  1995].  According  to  the  National 
Weather  Service  Weekly  Climate  Bulletin  (No.  92/37),  the  nation  as  a  whole  experienced 
its  third  coldest  and  third  wettest  summer  in  1992  since  records  began  in  1895.  1992  was 
the  nation’s  coolest  summer  in  77  years,  with  temperatures  below  normal  across  the 
eastern  two-thirds  of  the  nation.  The  South  experienced  the  second  coolest  summer  on 
record,  and  the  southern  Plains  and  Southeast  were  inundated  by  moderate  to  heavy  rains 
during  this  season. 

In  contrast  to  the  cooler  temperatures  and  above  normal  precipitation  during  1992, 
June  and  July  1993  was  dominated  by  a  persistent  circulation  pattern,  preventing  the 
eastward  progression  of  weather  systems  which  would  have  brought  rain  and  cooler 
weather  to  relieve  the  Southeast  of  drought  conditions  and  record  breaking  high 
temperatures.  The  Southeast  experienced  the  second  warmest  summer  in  99  years  of 
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record  [NWS  Weekly  Climate  Bulletin  93/39],  The  difference  in  ozone  levels  between 
1992  and  1993  is  explained  by  the  different  weather  patterns  that  dominated  each  summer. 
Ozone  is  positively  correlated  with  high  temperatures  and  negatively  correlated  with 
relative  humidity;  the  lower  ozone  in  1992  is  explained  by  the  relatively  lower 
temperatures  and  higher  relative  humidity,  as  well  as  the  greater  frequency  of  frontal 
passages,  which  tend  to  reduce  the  occurrence  of  high  pressure  systems  that  allow 
anthropogenic  pollutants  to  accumulate.  In  a  study  by  Vukovich  et  al.  [1977],  it  was 
found  that  high  ozone  in  the  summer  months  is  associated  with  high  pressure  systems,  and 
that  the  largest  concentrations  of  ozone  were  found  on  the  back  side  of  a  moving  high 
pressure  system  where  air  parcels  had  the  largest  residence  time.  A  modeling  study  for 
western  Australia  by  Hurley  and  Manins  [1995]  revealed  that  practically  all  high  ozone 
days  were  associated  with  recirculation  of  ozone  or  its  precursors. 

A  multiple  input  regression  with  ARMA(1,1)  errors  was  next  applied  to  the  1993 
data  to  compare  with  the  1992  regression  (Appendix  C).  The  regression  equation  came 
out  to  be 

Y,  =  1 .629T  -  0.539RH.  +  0. 1 59CO  +  0.29SO2,  +  1 .025NOy,  +  Z,  (7) 

where  Zt  =  0.949Zt.i  +  et  -  0.542et.i 

Figure  4a  shows  the  observed  ozone  overlaid  with  ozone  predicted  by  equation  7.  Figure 
4b  is  a  plot  of  residuals  with  predicted  ozone,  and  the  pattern  of  scatter  is  representative 
of  white  noise,  indicating  that  the  model  is  a  good  fit.  The  difference  between  the  1992 
regression  and  this  one  is  the  addition  of  temperature  and  carbon  monoxide,  and  the 
absence  of  barometric  pressure  in  the  1993  model.  This  suggests  a  greater  influence  of 
anthropogenic  emissions  on  ozone.  The  t-statistic  for  the  parameter  estimates  of  SO2  and 
NOy  shows  that  they  are  not  very  significant  in  the  model.  While  statistically  it  may  make 
sense  to  remove  them,  NOy  is  a  direct  precursor  to  photochemical  ozone  formation,  and  in 
an  atmospheric  chemistry  perspective  it  makes  more  sense  to  leave  it  in  the  model.  In  a 
separate  regression,  NOy  was  removed  and  the  t-statistic  of  CO  and  SO2  increased 
considerably  (Appendix  C,  Section  III). 
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Figure  4a.  Multiple  regression  of  Yorkville  1993  ozone  data  with  ARMA  errors 
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Figure  4b.  Predicted  ozone  versus  residuals  for  1993  regression 
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Summary 

A  multiple  regression  using  meteorological  and  trace  gas  data  as  input  parameters 
and  ARMA  time  series  errors  was  applied  to  model  daily  average  ozone  at  the  rural  site  of 
Yorkville,  Georgia,  using  data  collected  from  June  1  through  August  3 1  1992.  This 
statistical  model  was  then  applied  to  the  same  site  for  the  1993  summer.  It  was  found 
that,  while  a  good  regression  model  was  fitted  for  the  1992  ozone  data,  the  same  model 
tended  to  underpredict  observed  ozone  concentrations  for  the  1993  summer.  Ozone  was 
found  to  be  statistically  higher  in  1993  than  in  1992  (with  a  1993  average  of  63.5  ±  16.6 
ppb  and  a  1992  average  of  50.25  ±  14.8  ppb),  and  this  appears  to  be  explained  by  the 
synoptic  meteorology  characterizing  the  1993  summer,  which  was  conducive  to  high 
ozone  formation.  June  and  July  1993  was  dominated  by  a  persistent  circulation  pattern, 
preventing  the  eastward  progression  of  weather  systems  which  would  have  brought  rain 
and  cooler  weather  to  relieve  the  Southeast  of  record  breaking  high  temperatures  and 
drought  conditions.  In  contrast  to  the  stagnant  summer  conditions  of  the  1993  summer,  in 
1992  the  South  experienced  the  second  coolest  summer  on  record,  and  the  southern  Plains 
and  Southeast  were  inundated  by  moderate  to  heavy  rains  during  this  season.  In 
conclusion,  a  multiple  regression  model  with  time  series  errors  is  not  adequate  to  predict 
ozone  concentrations,  because  other  factors  affecting  ozone  levels,  such  as  synoptic 
meteorology  and  atmospheric  chemistry,  also  needs  to  be  considered. 
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APPENDICES 


Appendix  A.  Multiple  Regression  of  1992  Ozone  Data  with  ARMA  errors 


1.  Program 

goptions 
gunit=pct 
cback=white 
htitle=4 
htext=2 
ftext=swi3sb 
colors= (black) 
rotate=landscape; 

data  a; 

infile  'york92.dat'; 

input  Date  mmddyyS .  temp  relhum  pbaro  solrad  ozone  maxo3  CO  S02  NO  NOy; 
format  Date  date.; 

proc  arima  data=a; 

identify  var=02one  cros s cor r= (relhum  pbaro  S02  NOy)  noprint; 

estimate  p=l  q=l  input= (relhim  pbaro  S02  NOy)  noconstant  plot  method=ml; 

forecast  lead=0  id=Date  out=folder; 

proc  gplot  data=folder; 
title  'Yorkville  1992'; 
plot  (forecast  ozone) *Date/overlay; 
symboll  line=2  ; 
symbol2  line=l; 

proc  print  data=folder; 

run; 


11.  Statistical  Output  of  Yorkville  1992  Regression 


ARIMA  Procedure 
Maximum  Likelihood  Estimation 


Approx . 


Parcimeter 

Estimate 

Std  Error 

T  Ratio 

Lag 

Variable 

Shift 

MAI,  1 

0.66201 

0.16512 

4.01 

1 

OZONE 

0 

ARl,  1 

0.87511 

0.10515 

8.32 

1 

OZONE 

0 

NUMl 

-0.75793 

0.08486 

-8.93 

0 

RELHUM 

0 

NUM2 

0.13672 

0.0097887 

13.97 

0 

PBARO 

0 

NUM3 

-0.57661 

0.17002 

-3.39 

0 

S02 

0 

NUM4 

1.05332 

0.18372 

5.73 

0 

NOY 

0 

Variance 

Estimate  = 

79. 

4084332 

Std  Error 

Estimate  = 

8.91114096 

AIC 

= 

641 

.013713 

SBC 

= 

655 

.877734 

Number  of 

Residuals= 

88 
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Correlations  of  the  Estimates 


OZONE 

OZONE 

RELHUM 

PBARO 

S02 

Variable 

Parameter 

MA1,1 

ARl,  1 

NUMl 

NUM2 

NUM3 

OZONE 

MAI,  1 

1.000 

0.867 

0.073 

-0.109 

-0.042 

OZONE 

ARl,  1 

0.867 

1.000 

0.021 

-0.062 

-0.078 

RELHUM 

NUMl 

0.073 

0.021 

1.000 

-0.921 

0.017 

PBARO 

NXM2 

-0.109 

-0.062 

-0.921 

1.000 

-0.022 

S02 

NUM3 

-0.042 

-0.078 

0.017 
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1.000 

NOY 

NUM4 
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0.147 

0.136 
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NUMl 
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NUM2 

-0.288 

S02 

NUM3 

-0.515 

NOY 

NUM4 

1.000 

Autocorrelation  Check  of  Residuals 

To 
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Square 

DF 

Prob 

6 
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4 
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12 

11.92 

10 
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-0.035  0.109  -0.073 

0.004 
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18 

25.16 

16 
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-0.123 

24 

32.81 

22 

0.065 

-0.173 

-0.143  -0.056  -0.014 

-0.059 

-0.082 

Autocorrelation  Plot  of  Residuals 
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Correlation  -1987654321 

01234567891 

0 
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1.00000  1 

1  ******************** 

1 
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0.05375  1 

1  * 

2 
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1 

3 
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* 

4 
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1 

5 
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6 
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1 

7 
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8 
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★ 

9 
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1  *  *  ^ 

10 

-5.437432 

-0.07348  1 

♦ 

11 

0.319807 

0.00432  1 

12 
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I  *  *  *  *  ^ 

13 

-13.190454 

-0.17826  1 

14 

4.731964 
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1  * 

15 

-12.325864 
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^  * 

16 

-5.342739 
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★ 

17 

13.670529 
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18 
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** 

19 

-12.796916 

-0.17294  1 

^  ★ 

20 

-10.598759 

-0.14324  1 

^  ★ 

21 

-4.175482 

-0.05643  1 

★ 

22 

-1.061787 

-0.01435  1 

23 

-4.365145 

-0.05899  1 

* 

marks  two  standard  errors 
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Model  for  variable  OZONE 


No  mean  term  in  this  model. 


Autoregressive  Factors 
Factor  1:  1  -  0.87511  B**(l) 

Moving  Average  Factors 
Factor  1:  1  -  0.66201  B**(l) 


Input  Number  1  is  RELHUM. 
Overall  Regression  Factor 

Input  Number  2  is  PBARO. 
Overall  Regression  Factor 

Input  Number  3  is  S02. 
Overall  Regression  Factor 

Input  Number  4  is  NOY. 
Overall  Regression  Factor 


=  -0.75793 

=  0.136722 

=  -0.57661 

=  1.053324 


I/I.  Output  of  observed  and  predicted  ozone,  standard  deviation,  upper  and  lower  95% 
confidence  intervals  and  residuals. 

Yorkville  1992  Output 


OBS 

DATE 

OZONE 

FORECAST 

STD 

L95 

U95 

RESIDUAL 

1 

01JUN92 

46.46 

52.5920 

9.73688 

33.5080 

71.6759 

-6.1320 

2 

02JUN92 

81.46 

67.0483 

9.22284 

48.9719 

85.1247 

14.4117 

3 

03JUN92 

35.76 

29.3509 

9.03997 

11.6329 

47.0689 

6.4091 

4 

04JUN92 

57.04 

52.1065 

8.96623 

34.5330 

69.6800 

4.9335 

5 

05JUN92 

55.34 

58.4955 

8.93503 

40.9832 

76.0079 

-3.1555 

6 

06JUN92 

65.40 

70.3747 

8.92156 

52.8888 

87.8607 

-4.9747 

7 

07Jtm92 

49.70 

53.3187 

8.91570 

35.8443 

70.7932 

-3.6187 

8 

08JUN92 

18.81 

28.9204 

8.91314 

11.4510 

46.3898 

-10.1104 

9 

09JUN92 

24.18 

29.8179 

8.91202 

12.3507 

47.2851 

-5.6379 

10 

10JUN92 

48.88 

45.5842 

8.91152 

28.1179 

63.0505 

3.2958 

11 

11JUN92 

57.26 

42.5497 

8.91131 

25.0839 

60.0156 

14.7103 

12 

12JUN92 

44.92 

46.8236 

8.91121 

29.3579 

64.2892 

-1.9036 

13 

13JUN92 

24.78 

40.6999 

8.91117 

23.2344 

58.1655 

-15.9199 

14 

14JUN92 

40.47 

40.4147 

8.91117 

22.9491 

57.8802 

0.0553 

15 

15JUN92 

52.13 

47.3188 

8.91117 

29.8532 

64.7844 

4.8112 

16 

16JUN92 

81.83 

59.2664 

8.91117 

41.8008 

76.7320 

22.5636 

17 

17JUN92 

56.92 

57.7109 

8.91117 

40.2454 

75.1765 

-0.7909 

18 

18JUN92 

41.67 

48.3452 

8.91117 

30.8796 

65.8108 

-6.6752 

19 

19JUN92 

65.62 

61.0159 

8.91117 

43.5503 

78.4815 

4.6041 

20 

20JUN92 

67.74 

64.2845 

8.91117 

46.8189 

81.7501 

3.4555 

21 

21JUN92 

52.41 

51.0400 

8.91117 

33.5744 

68.5056 

1.3700 

22 

22JUN92 

68.88 

72.9288 

8.91117 

55.4632 

90.3944 

-4.0488 

23 

23JUN92 

81.63 

66.5668 

8.91117 

49.1013 

84.0324 

15.0632 

24 

24JUN92 

64.81 

74.2631 

8.91117 

56.7975 

91.7287 

-9.4531 

25 

25JUN92 

64.44 

63.6685 

8.91117 

46.2029 

81.1341 

0.7715 

26 

26JUN92 

58.90 

50.8632 

8.91117 

33.3977 

68.3288 

8.0368 

27 

27JUN92 

57.88 

53.0958 

8.91117 

35.6302 

70.5614 

4.7842 

28 

28JUN92 

75.90 

57.2029 

8.91117 

39.7374 

74.6685 

18.6971 

29 

29JUN92 

76.08 

69.0518 

8.91117 

51.5862 

86.5174 

7.0282 

30 

30JUN92 

46.46 

46.4082 

8.91117 

28.9426 

63.8738 

0.0518 

31 

01JUL92 

33.26 

41.3676 

8.91117 

23.9021 

58.8332 

-8.1076 

32 

02JUL92 

, 

47.6725 

8.91117 

30.2070 

65.1381 

33 

03JUL92 

37.81 

43.3168 

9.11125 

25.4591 

61.1745 

-5.5068 

34 

04JUL92 

56.01 

55.9081 

8.99557 

38.2771 

73.5391 

0.1019 

35 

05JUL92 

41.43 

61.1407 

8.94755 

43.6039 

78.6776 

-19.7107 
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36 

06JUL92 

41.42 

42.7796 

8.92699 

25.2831 

60.2762 

-1.3596 

37 

07JUL92 

64.87 

54.6953 

8.91806 

37.2163 

72.1744 

10.1747 

38 

08JUL92 

, 

73.5462 

8.91417 

56.0748 

91.0177 

, 

39 

09JUL92 

48.70 

59.9744 

9.11350 

42.1123 

77.8366 

-11.2744 

40 

10JUL92 

57.44 

56.0101 

8.99649 

38.3773 

73.6428 

1.4299 

41 

11JUL92 

49.69 

63.0395 

8.94794 

45.5019 

80.5771 

-13.3495 

42 

12JUL92 

53.78 

64.0271 

8.92715 

46.5302 

81.5240 

-10.2471 

43 

13JUL92 

53.47 

45.2427 

8.91814 

27.7635 

62.7219 

8.2273 

44 

14JUL92 

40.68 

50.7310 

8.91420 

33.2595 

68.2025 

-10.0510 

45 

15JUL92 

34.75 

30.8658 

8.91248 

13.3977 

48.3340 

3.8842 

46 

16JUL92 

35.67 

41.5735 

8.91173 

24.1068 

59.0402 

-5.9035 

47 

17JUL92 

31.15 

40.0879 

8.91140 

22.6219 

57.5539 

-8.9379 

48 

18JUL92 

39.78 

33.1792 

8.91125 

15.7134 

50.6449 

6.6008 

49 

19JUL92 

56.23 

70.7538 

8.91119 

53.2882 

88.2194 

-14.5238 

50 

20JUL92 

74.16 

51 

21JUL92 

49.04 

51.3692 

9.11125 

33.5115 

69.2269 

-2.3292 

52 

22JUL92 

35.71 

30.8028 

8.99557 

13.1718 

48.4338 

4.9072 

53 

23JUL92 

31.52 

39.4495 

8.94755 

21.9126 

56.9864 

-7.9295 

54 

24JUL92 

33.04 

38.3068 

8.92699 

20.8102 

55.8034 

-5.2668 

55 

25JUL92 

40.65 

50.5007 

8.91806 

33.0217 

67.9798 

-9.8507 

56 

26JUL92 

38.93 

45.7602 

8.91417 

28.2887 

63.2316 

-6.8302 

57 

27JUL92 

46.36 

36.0857 

8.91247 

18.6176 

53.5538 

10.2743 

58 

28JUL92 

48.56 

47.5803 

8.91172 

30.1136 

65.0469 

0.9797 

59 

29JUL92 

57.61 

53.0988 

8.91140 

35.6328 

70.5648 

4.5112 

60 

30JUL92 

47.28 

43.1727 

8.91125 

25.7069 

60.6384 

4.1073 

61 

31JUL92 

37.46 

42.0449 

8.91119 

24.5793 

59.5105 

-4.5849 

62 

01AUG92 

51.31 

52.3051 

8.91116 

34.8395 

69.7706 

-0.9951 

63 

02AUG92 

52.54 

61.3151 

8.91116 

43.8495 

78.7807 

-8.7751 

64 

03AUG92 

60.02 

58.2078 

8.91116 

40.7422 

75.6733 

1.8122 

65 

04AUG92 

53.06 

51.0435 

8.91116 

33.5779 

68.5090 

2.0165 

66 

05AUG92 

61.82 

58.7871 

8.91116 

41.3215 

76.2527 

3.0329 

67 

06AUG92 

85.78 

77.6196 

8.91116 

60.1541 

95.0852 

8.1604 

68 

07AUG92 

44.76 

39.5477 

8.91116 

22.0822 

57.0133 

5.2123 

69 

08AUG92 

55.21 

42.4226 

8.91116 

24.9571 

59.8882 

12.7874 

70 

09AUG92 

49.92 

52.9368 

8.91116 

35.4713 

70.4024 

-3.0168 

71 

10AUG92 

54.35 

56.9933 

8.91116 

39.5277 

74.4588 

-2.6433 

72 

11AUG92 

48.51 

58.3014 

8.91116 

40.8358 

75.7669 

-9.7914 

73 

12AUG92 

48.66 

36.0146 

8.91116 

18.5491 

53.4802 

12.6454 

74 

13AUG92 

35.30 

27.4822 

8.91116 

10.0167 

44.9478 

7.8178 

75 

14AUG92 

39.30 

42.6036 

8.91116 

25.1381 

60.0692 

-3.3036 

76 

15AUG92 

56.81 

41.3615 

8.91116 

23.8960 

58.8271 

15.4485 

77 

16AUG92 

48.20 

41.6379 

8.91116 

24.1723 

59.1034 

6.5621 

78 

17AUG92 

39.46 

44.8872 

8.91116 

27.4217 

62.3528 

-5.4272 

79 

18AUG92 

51.76 

48.6345 

8.91116 

31.1689 

66.1000 

3.1255 

80 

19AUG92 

69.52 

56.4268 

8.91116 

38.9613 

73.8924 

13.0932 

81 

20AUG92 

73.20 

62.6012 

8.91116 

45.1357 

80.0668 

10.5988 

82 

21AUG92 

25.76 

42.0489 

8.91116 

24.5833 

59.5145 

-16.2889 

83 

22AUG92 

33.71 

35.2559 

8.91116 

17.7904 

52.7215 

-1.5459 

84 

23AUG92 

34.42 

41.4708 

8.91116 

24.0052 

58.9363 

-7.0508 

85 

24AUG92 

55.33 

50.5741 

8.91116 

33.1086 

68.0397 

4.7559 

86 

25AUG92 

62.43 

54.0708 

8.91116 

36.6053 

71.5364 

8.3592 

87 

26AUG92 

26.35 

42.2938 

8.91116 

24.8283 

59.7594 

-15.9438 

88 

27AUG92 

23.83 

30.6095 

8.91116 

13.1440 

48.0751 

-6.7795 

89 

28AUG92 

25.81 

35.2953 

8.91116 

17.8297 

52.7608 

-9.4853 

90 

29AUG92 

41.32 

43.9097 

8.91116 

26.4442 

61.3753 

-2.5897 

91 

30AUG92 

57.32 

50.4844 

8.91116 

33.0188 

67.9499 

6.8356 
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Appendix  B.  Application  of  1992  model  for  1993  Ozone  Data 

1.  Program 

goptions 
giinit=pct 
cback=white 
htitle=4 
htext=2 
ftext=swissb 
colors= (black) 
rotate=landscape; 

data  a; 

infile  'york93.dat'; 

input  Date  nimddyyS.  temp  relhum  pbaro  solrad  ozone  maxo3  CO  S02  NO  NOy 
format  Date  date.; 

proc  arima  data=a; 

identify  var=ozone  crosscorr= (relhum  pbaro  S02  NOy)  noprint; 
estimate  p=l  q=l  input= (relhum  pbaro  S02  NOy) 

ar=0. 87511  ma=0. 66201  ini tval= (-0 . 75793  relhum 

0.136722  pbaro  -0.57661  S02  1.053324  NOy  )  noconstant  noest 

method=ml; 

forecast  lead=0  id=Date  out=folder; 
proc  print  data=folder; 
proc  gplot; 

title  'Yorkville  1993'; 

plot  (forecast  ozone) *Date/overlay; 

symbol 1  line=2; 

3ymbol2  line=l; 

proc  univariate  data=folder; 
var  residual; 

run; 


11.  Statistical  output  of 1992  regression  model  on  1993  ozone  data 


ARIMA  Procedure 


Maximum  Likelihood  Estimation 


Parameter 

MAI,  1 

ARl,  1 

NUMl 

NUM2 

NUM3 

NUM4 


Estimate 

0.66201 

0.87511 

-0.75793 

0.13672 

-0.57661 

1.05332 


Approx . 

Std  Error  T  Ratio 
0 
0 
0 
0 
0 
0 


Lag  Variable  Shift 

1  OZONE  0 

1  OZONE  0 

0  RELHUM  0 

0  PBARO  0 

0  S02  0 

0  NOY  0 


Variance  Estimate  = 
Std  Error  Estimate  = 
AIC 
SBC 

Number  of  Residuals= 


195.664503 

13.9880128 

736.372255 

751.371113 
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Autocorrelation  Check  of  Residuals 


To 

Chi 

Autocorrelations 

Lag 

Square 

DF 

Prob 

6 

21.41 

4 

0.000 

0.215 

0.088 

0.095 

0.230 

0.288 

0.154 

12 

30.80 

10 

0.001 

0.114 

0.086 

0.156 

0.055 

0.202 

0.047 

18 

35.23 

16 

0.004 

0.070 

0.033 

0.115 

0.108 

-0.086 

-0.041 

24 

39.82 

22 

0.011 

0.063 

0.118 

-0.025 

0.048 

0.001 

0.130 

HI.  Output  of  observed  and  predated  ozone,  standard  deviation,  upper  and  lower  95% 
confidence  intervals  and  residuals. 


OBS 

DATE 

OZONE 

FORECAST 

STD 

L95 

U95 

RESIDUAL 

1 

01Jra93 

55.04 

2 

02JUN93 

68.24 

. 

3 

03JUN93 

71.48 

58.6499 

15.2842 

28.6934 

88.6064 

12.8301 

4 

04JUN93 

36.20 

50.8546 

14.4773 

22.4796 

79.2296 

-14.6546 

5 

05JUN93 

44.04 

35.6236 

14.1902 

7.8113 

63.4360 

8.4164 

6 

06JUN93 

57.28 

61.3149 

14.0745 

33.7294 

88.9003 

-4.0349 

7 

07JUN93 

62.52 

52.7599 

14.0255 

25.2705 

80.2494 

9.7601 

8 

08JUN93 

50.84 

56.5240 

14.0044 

29.0759 

83.9721 

-5.6840 

9 

09JUN93 

42.20 

56.4251 

13.9952 

28.9950 

83.8551 

-14.2251 

10 

10JUN93 

53.80 

55.0555 

13.9911 

27.6334 

82.4777 

-1.2555 

11 

11JUN93 

52.28 

55.9915 

13.9894 

28.5728 

83.4102 

-3.7115 

12 

12JUN93 

84.46 

40.9817 

13.9886 

13.5645 

68.3989 

43.4783 

13 

13JUN93 

71.88 

54.9182 

13.9883 

27.5017 

82.3347 

16.9618 

14 

14JUN93 

57.32 

46.9154 

13.9881 

19.4992 

74.3316 

10.4046 

15 

15JUN93 

65.12 

60.6445 

13.9880 

33.2285 

88.0605 

4.4755 

16 

16JUN93 

81.56 

61.4937 

13.9880 

34.0777 

88.9097 

20.0663 

17 

17JUN93 

76.84 

70.7463 

13.9880 

43.3303 

98.1623 

6.0937 

18 

18JUN93 

77.20 

71.1590 

13.9880 

43.7430 

98.575 

6.0410 

19 

19JUN93 

73.36 

73.1400 

13.9880 

45.7240 

100.556 

0.2200 

20 

20JtJN93 

58.52 

62.1667 

13.9880 

34.7507 

89.583 

-3.6467 

21 

21JUN93 

36.36 

37.9696 

13.9880 

10.5536 

65.386 

-1.6096 

22 

22JUN93 

32.67 

46.4388 

13.9880 

19.0228 

73.855 

-13.7688 

23 

23JUN93 

93.76 

63.6147 

13.9880 

36.1987 

91.031 

30.1453 

24 

24JUN93 

65.64 

59.2929 

13.9880 

31.8769 

86.709 

6.3471 

25 

25JUN93 

41.40 

37.8363 

13.9880 

10.4203 

65.252 

3.5637 

26 

26JUN93 

44.52 

54.0530 

13.9880 

26.6370 

81.469 

-9.5330 

27 

27JUN93 

62.28 

56.6642 

13.9880 

29.2482 

84.080 

5.6158 

28 

28JUN93 

44.16 

51.5753 

13.9880 

24.1593 

78.991 

-7.4153 

29 

29JUN93 

42.72 

41.5124 

13.9880 

14.0964 

68.928 

1.2076 

30 

30JUN93 

44.08 

49.2412 

13.9880 

21.8252 

76.657 

-5.1612 

31 

01JUL93 

46.92 

49.2137 

13.9880 

21.7977 

76.630 

-2.2937 

32 

02JUL93 

56.76 

51.2595 

13.9880 

23.8435 

78.676 

5.5005 

33 

03JUL93 

47.28 

59.5192 

13.9880 

32.1032 

86.935 

-12.2392 

34 

04JUL93 

54.60 

59.2284 

13.9880 

31.8124 

86.644 

-4.6284 

35 

05JUL93 

40.04 

45.9888 

13.9880 

18.5728 

73.4048 

-5.9488 

36 

06JUL93 

37.68 

58.4776 

13.9880 

31.0616 

85.8936 

-20.7976 

37 

07JUL93 

55.20 

55.1551 

13.9880 

27.7391 

82.5711 

0.0449 

38 

08JUL93 

73.80 

50.5884 

13.9880 

23.1724 

78.0044 

23.2116 

39 

09JUL93 

69.16 

57.7131 

13.9880 

30.2971 

85.1291 

11.4469 

40 

10JUL93 

55.52 

59.4458 

13.9880 

32.0298 

86.8618 

-3.9258 

41 

11JUL93 

41.24 

60.3920 

13.9880 

32.9760 

87.8080 

-19.1520 

42 

12JUL93 

50.96 

51.1649 

13.9880 

23.7489 

78.5809 

-0.2049 

43 

13JUL93 

63.00 

36.1074 

13.9880 

8.6913 

63.5234 

26.8926 

44 

14JUL93 

41.88 

55.9865 

13.9880 

28.5705 

83.4025 

-14.1065 

45 

15JUL93 

48.56 

56.0475 

13.9880 

28.6315 

83.4635 

-7.4875 

46 

16JUL93 

45.44 

37.0965 

13.9880 

9.6805 

64.5125 

8.3435 

47 

17JUL93 

69.08 

50.8909 

13.9880 

23.4749 

78.3069 

18.1891 

48 

18JUL93 

76.84 

61.4443 

13.9880 

34.0283 

88.8603 

15.3957 

49 

19JUL93 

52.32 

64.2552 

13.9880 

36.8392 

91.6712 

-11.9352 

50 

20JUL93 

69.60 

60.5386 

13.9880 

33.1226 

87.9546 

9.0614 

51 

21JUL93 

88.04 

70.2391 

13.9880 

42.8231 

97.6551 

17.8009 
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52 

22JUL93 

95.12 

68.4536 

13.9880 

41.0376 

95.870 

26.6664 

53 

23JUL93 

91.60 

70.5522 

13.9880 

43.1362 

97.968 

21.0478 

54 

24JUL93 

81.40 

68.8673 

13.9880 

41.4513 

96.283 

12.5327 

55 

25JUL93 

57.20 

63.5778 

13.9880 

36.1618 

90.994 

-6.3778 

56 

26JUL93 

74.96 

65.5030 

13.9880 

38.0870 

92.919 

9.4570 

57 

27JUL93 

74.67 

69.1706 

13.9880 

41.7546 

96.587 

5.4994 

58 

28JUL93 

89.72 

71.7687 

13.9880 

44.3527 

99.185 

17.9513 

59 

29JUL93 

93.48 

74.9488 

13.9880 

47.5328 

102.365 

18.5312 

60 

30JUL93 

64.12 

84.6861 

13.9880 

57.2701 

112.102 

-20.5661 

61 

31JUL93 

71.16 

82.6447 

13.9880 

55.2287 

110.061 

-11.4847 

62 

01AUG93 

82.04 

69.8863 

13.9880 

42.4703 

97.302 

12.1537 

63 

02AUG93 

68.44 

60.8106 

13.9880 

33.3946 

88.227 

7.6294 

64 

03AUG93 

48.87 

45.3704 

13.9880 

17.9544 

72.786 

3.4996 

65 

04AUG93 

50.40 

42.6967 

13.9880 

15.2807 

70.113 

7.7033 

66 

05AUG93 

55.88 

56.3689 

13.9880 

28.9529 

83.785 

-0.4889 

67 

06AUG93 

39.00 

34.4386 

13.9880 

7.0226 

61.855 

4.5614 

68 

07AUG93 

35.84 

35.0836 

13.9880 

7.6676 

62.500 

0.7564 

69 

08AUG93 

51.12 

49.8953 

13.9880 

22.4793 

77.311 

1.2247 

70 

09AUG93 

69.64 

54.6847 

13.9880 

27.2687 

82.101 

14.9553 

71 

10AUG93 

71.88 

63.4367 

13.9880 

36.0207 

90.853 

8.4433 

72 

11AUG93 

76.55 

77.8721 

13.9880 

50.4561 

105.288 

-1.3221 

73 

12AUG93 

64.08 

65.5517 

13.9880 

38.1357 

92.968 

-1.4717 

74 

13AUG93 

52.08 

45.2724 

13.9880 

17.8564 

72.688 

6.8076 

75 

14AUG93 

55.24 

47.3098 

13.9880 

19.8938 

74.726 

7.9302 

76 

15AUG93 

63.16 

51.8960 

13.9880 

24.4800 

79.312 

11.2640 

77 

16AUG93 

75.00 

55.8594 

13.9880 

28.4434 

83.275 

19.1406 

78 

17AUG93 

77.52 

63.6965 

13.9880 

36.2805 

91.113 

13.8235 

79 

18AUG93 

76.24 

59.8057 

13.9880 

32.3897 

87.222 

16.4343 

80 

19AUG93 

101.64 

67.5072 

13.9880 

40.0912 

94.923 

34.1328 

81 

20AUG93 

75.96 

75.0659 

13.9880 

47.6499 

102.482 

0.8941 

82 

21AUG93 

74.12 

64.4080 

13.9880 

36.9920 

91.824 

9.7120 

83 

22AUG93 

70.60 

69.6054 

13.9880 

42.1894 

97.021 

0.9946 

84 

23AUG93 

93.84 

67.4777 

13.9880 

40.0617 

94.894 

26.3623 

85 

24AUG93 

57.42 

68.2805 

13.9880 

40.8645 

95.697 

-10.8605 

86 

25AUG93 

77.28 

61.5506 

13.9880 

34.1346 

88.967 

15.7294 

87 

26AUG93 

101.68 

75.2405 

13.9880 

47.8245 

102.656 

26.4395 

88 

27AUG93 

67.52 

76.1929 

13.9880 

48.7769 

103.609 

-8.6729 

89 

28AUG93 

73.64 

68.8648 

13.9880 

41.4488 

96.281 

4.7752 

90 

29AUG93 

70.88 

67.8343 

13.9880 

40.4183 

95.250 

3.0457 

91 

30AUG93 

60.84 

69.7646 

13.9880 

42.3486 

97.181 

-8.9246 

92 

31AUG93 

75.96 

67.8189 

13.9880 

40.4029 

95.235 

8.1411 

IV.  Descriptive  statistics  of  the  residuals 


Univariate  Procedure 


Variable=RESIDUAL  Residual:  Actual-Forecast 


Moments 


N 

90 

Mean 

4.91295 

Std  Dev 

12.67913 

Skewness 

0.332227 

USS 

16480 

CV 

258.0757 

T:Mean=0 

3.675989 

Num  "=  0 

90 

M{Sign) 

13 

Sgn  Rank 

810.5 

Sum  Wgts 

90 

Sum 

442.1655 

Variance 

160.7603 

Kurtosis 

0.156777 

CSS 

14307.67 

Std  Mean 

1.336497 

Pr>|T| 

0.0004 

Num  >  0 

58 

Pr>=|M| 

0.0080 

Pr>= 1 S 1 

0.0008 
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Appendix  C.  Multiple  Regression  on  1993  ozone  data  with  ARMA  errors 


1.  Program 


goptions 
gunit=pct 
cback=white 
htitle=4 
htext=2 
ftext=swissb 
colors= (black) 
rotate=landscape; 

data  a; 

infile  'york93.dat'; 

input  Date  mmddyyS .  temp  relhum  pbaro  solrad  ozone  maxo3  CO  S02  NO  NOy; 
format  Date  date.; 

proc  arima; 

identify  var=o2one  cro3scorr= (temp  relhum  pbaro  solrad  CO  S02  NOy) 
noprint; 

estimate  p=l  q=l  input=(temp  relhim  CO  S02  NOy) 
noconstant  plot  method=ml; 
forecast  lead=0  id=Date  out=folder; 

proc  arima  data=a; 

identify  var=ozone  crosscorr= (temp  relhum  pbaro  solrad  CO  S02  NOy) 
noprint ; 

estimate  p=l  q=l  input=(temp  relhum  CO  S02) 
noconstant  plot  method=ml; 
forecast  lead=0  id=Date  out=folder2; 

proc  gplot  data=folder; 

title  'Yorkville  1993'; 

plot  (forecast  ozone) *Date/overlay; 

symbol 1  line=2; 

symbol2  line=l; 

proc  print  data=folder; 

run; 


11.  Statistical  output  of 1993  regression  (NOy  is  included  in  regression) 


ARIMA  Procedure 


Maximum  Likelihood  Estimation 


Parameter 

Estimate 

Approx . 
Std  Error 

T  Ratio 

Lag 

Variab] 

MAI,  1 

0.54244 

0.11502 

4.72 

1 

OZONE 

ARl,  1 

0.94866 

0.04309 

22.01 

1 

OZONE 

NUMl 

1.62859 

0.21083 

7.72 

0 

TEMP 

NUM2 

-0.53881 

0.07426 

-7.26 

0 

RELHUM 

NUM3 

0.15885 

0.02271 

7.00 

0 

CO 

NUM4 

0.28976 

0.14814 

1.96 

0 

S02 

NUM5 

1.02496 

0.59466 

1.72 

0 

NOY 
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Variance  Estimate  = 
Std  Error  Estimate  = 
AIC 
SBC 

Number  of  Residuals= 


51.0245468 

7.14314684 

617.234365 

634.733033 

90 


Correlations  of  the  Estimates 


OZONE 

OZONE 

TEMP 

RELHUM 

CO 

Variable 

Parameter 

MAI,  1 

ARl,  1 

NUMl 

NUM2 

NUM3 

OZONE 

MAI,  1 

1.000 

0.582 

0.012 

0.080 

-0.043 

OZONE 

ARl,  1 

0.582 

1.000 

0.020 

-0.107 

0.142 

TEMP 

NUMl 

0.012 

0.020 

1.000 

-0.057 

-0.264 

RELHUM 

NUM2 

0.080 

-0.107 

0.057 

1.000 

-0.643 

CO 

NUM3 

-0.043 

0.142 

0.264 

-0.643 

1.000 

S02 

NUM4 

-0.087 

0.090 

0.102 

-0.267 

0.479 

NOY 

NUM5 

0.092 

-0.120 

0.014 

0.520 

-0.761 

S02 

NOY 

Variable 

Parameter  NUM4 

NUM5 

OZONE 

MAI,  1 

-0.087 

0.092 

OZONE 

ARl,  1 

0.090 

- 

0.120 

TEMP 

NUMl 

-0.102 

- 

0.014 

RELHUM 

NUM2 

-0.267 

0.520 

CO 

NUM3 

0.479 

- 

0.761 

S02 

NUM4 

1.000 

- 

0.653 

NOY 

NUM5 

-0.653 

1.000 

Autocorrelation  Check  of  Residuals 
To  Chi  Autocorrelations 


Lag 

6 

Square 

2.50 

DF 

4 

Prob 

0.645 

-0.012 

-0.006 

0.103 

0.016 

0.121 

0.008 

12 

8.62 

10 

0.568 

-0.034 

0.002 

-0.225 

-0.063 

0.058 

0.025 

18 

16.28 

16 

0.434 

0.043 

-0.170 

-0.010 

0.183 

0.033 

0.060 

24 

20.43 

22 

0.556 

0.097 

0.098 

0.012 

-0.091 

0.078 

0.032 

Model  for  variable  OZONE 

No  mean  term  in  this  model. 

Autoregressive  Factors 
Factor  1:  1  -  0.94866  B**(l) 

Moving  Average  Factors 
Factor  1:  1  -  0.54244  B**(l) 


Input  Number  1  is  TEMP. 
Overall  Regression  Factor 

Input  Number  2  is  RELHUM. 
Overall  Regression  Factor 

Input  Number  3  is  CO. 
Overall  Regression  Factor 

Input  Number  4  is  S02 . 
Overall  Regression  Factor 

Input  Number  5  is  NOY. 
Overall  Regression  Factor 


=  1.628593 

=  -0.53881 

=  0.158854 

=  0.289762 

=  1.024958 
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HI.  Statistical  output  of  multiple  regression  on  the  1993  ozone  data  with  NOy  removed 


ARIMA  Procedure 


Maximum  Likelihood  Estimation 


Approx . 


Parameter 

Estimate 

Std  Error 

T  Ratio 

Lag 

Variable 

Shift 

MAI,  1 

0.53400 

0.11144 

4.79 

1 

OZONE 

0 

ARl,  1 

0.95798 

0.03840 

24.95 

1 

OZONE 

0 

NUMl 

1.61859 

0.23266 

6.96 

0 

TEMP 

0 

NUM2 

-0.60666 

0.06564 

-9.24 

0 

RELHUM 

0 

NUM3 

0.18906 

0.01480 

12.78 

0 

CO 

0 

NUM4 

0.45282 

0.11242 

4.03 

0 

S02 

0 

Variance  Estimate  =  52.0935906 


Std  Error  Estimate  =  7.21758897 
AIC  =  618.373558 
SBC  =  633.372416 
Number  of  Residuals=  90 


Autocorrelation  Check  of  Residuals 


To 

Chi 

Autocorrelations 

Lag 

Square 

DF 

Prob 

6 

3.63 

4 

0.458 

-0.017 

-0.008 

0.143 

-0.001 

0.130 

0.001 

12 

10.32 

10 

0.413 

-0.024 

0.034 

-0.243 

-0.038 

0.043 

-0.038 

18 

19.68 

16 

0.235 

0.067 

-0.176 

-0.028 

0.204 

-0.006 

0.082 

24 

23.52 

22 

0.373 

0.054 

0.123 

0.062 

-0.075 

0.052 

0.043 

Model  for  variable  OZONE 
No  mean  term  in  this  model. 


Autoregressive  Factors 
Factor  1:  1  -  0.95798  B**(l) 

Moving  Average  Factors 
Factor  1:  1  -  0.534  B**(l) 


Input  Number  1  is  TEMP. 

Overall  Regression  Factor  =  1.618591 


Input  Number  2  is  RELHUM. 
Overall  Regression  Factor 

Input  Number  3  is  CO. 
Overall  Regression  Factor 

Input  Number  4  is  S02. 
Overall  Regression  Factor 


-0.60666 

0.189061 

0.452822 
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IV.  Observed  and  predicted  ozone,  standard  deviation,  95%  confidence  intervals  and 
residuals  (from  1993  regression  that  includes  the  NOy  variable) 


OBS 

DATE 

OZONE 

FORECAST 

STD 

L95 

U95 

RESIDUAL 

1 

01JUN93 

55.04 

2 

02JUN93 

68.24 

3 

03JUN93 

71.48 

65.8940 

11.6266 

43.1063 

88.6817 

5.5860 

4 

04JUN93 

36.20 

55.6115 

7.7699 

40.3828 

70.8402 

-19.4115 

5 

05JUN93 

44.04 

41.2274 

7.3040 

26.9118 

55.5430 

2.8126 

6 

06JUN93 

57.28 

53.9778 

7.1888 

39.8880 

68.0676 

3.3022 

7 

07JUN93 

62.52 

70.5505 

7.1564 

56.5242 

84.5769 

-8.0305 

8 

08JUN93 

50.84 

56.6778 

7.1470 

42.6699 

70.6858 

-5.8378 

9 

09JUN93 

42.20 

49.0256 

7.1443 

35.0230 

63.0282 

-6.8256 

10 

10JUN93 

53.80 

52.0610 

7.1435 

38.0600 

66.0620 

1.7390 

11 

11JUN93 

52.28 

54.5785 

7.1432 

40.5780 

68.5791 

-2.2985 

12 

12JUN93 

84.46 

76.1472 

7.1431 

62.1469 

90.1475 

8.3128 

13 

13JUN93 

71.88 

63.9261 

7.1431 

49.9258 

77.9264 

7.9539 

14 

14JUN93 

57.32 

71.9026 

7.1431 

57.9022 

85.9029 

-14.5826 

15 

15JUN93 

65.12 

56.9703 

7.1431 

42.9700 

70.9706 

8.1497 

16 

16JUN93 

81.56 

73.7138 

7.1431 

59.7135 

87.7141 

7.8462 

17 

17JUN93 

76.84 

74.6258 

7.1431 

60.6255 

88.6261 

2.2142 

18 

18JUN93 

77.20 

73.4532 

7.14315 

59.4529 

87.454 

3.7468 

19 

19JUN93 

73.36 

69.2852 

7.14315 

55.2849 

83.285 

4.0748 

20 

20JUN93 

58.52 

70.6697 

7.14315 

56.6694 

84.670 

-12.1497 

21 

21JUN93 

36.36 

34.2456 

7.14315 

20.2453 

48.246 

2.1144 

22 

22JUN93 

32.67 

40.8658 

7.14315 

26.8655 

54.866 

-8.1958 

23 

23JUN93 

93.76 

86.4889 

7.14315 

72.4886 

100.489 

7.2711 

24 

24JUN93 

65.64 

72.3113 

7.14315 

58.3109 

86.312 

-6.6713 

25 

25JUN93 

41.40 

53.5279 

7.14315 

39.5276 

67.528 

-12.1279 

26 

26JUN93 

44.52 

43.6802 

7.14315 

29.6799 

57.681 

0.8398 

27 

27JUN93 

62.28 

60.3573 

7.14315 

46.3570 

74.358 

1.9227 

28 

28JUN93 

44.16 

43.3208 

7.14315 

29.3205 

57.321 

0.8392 

29 

29JUN93 

42.72 

43.6092 

7.14315 

29.6088 

57.609 

-0.8892 

30 

30JUN93 

44.08 

51.2911 

7.14315 

37.2908 

65.291 

-7.2111 

31 

01JUL93 

46.92 

48.9135 

7.14315 

34.9132 

62.914 

-1.9935 

32 

02JUL93 

56.76 

52.4943 

7.14315 

38.4940 

66.495 

4.2657 

33 

03JUL93 

47.28 

54.0724 

7.14315 

40.0721 

68.073 

-6.7924 

34 

04JUL93 

54.60 

60.7313 

7.14315 

46.7309 

74.732 

-6.1313 

35 

05JUL93 

40.04 

39.4462 

7.14315 

25.4459 

53.4465 

0.5938 

36 

06JUL93 

37.68 

51.8667 

7.14315 

37.8664 

65.8670 

-14.1867 

37 

07JUL93 

55.20 

58.7050 

7.14315 

44.7047 

72.7053 

-3.5050 

38 

08JUL93 

73.80 

73.9004 

7.14315 

59.9001 

87.9008 

-0.1004 

39 

09JUL93 

69.16 

61.3399 

7.14315 

47.3396 

75.3402 

7.8201 

40 

10JUL93 

55.52 

58.8754 

7.14315 

44.8751 

72.8758 

-3.3554 

41 

11JUL93 

41.24 

44.3816 

7.14315 

30.3813 

58.3819 

-3.1416 

42 

12JUL93 

50.96 

47.4887 

7.14315 

33.4884 

61.4890 

3.4713 

43 

13JUL93 

63.00 

61.0294 

7.14315 

47.0291 

75.0297 

1.9706 

44 

14JUL93 

41.88 

43.8170 

7.14315 

29.8167 

57.8173 

-1.9370 

45 

15JUL93 

48.56 

48.4649 

7.14315 

34.4646 

62.4652 

0.0951 

46 

16JUL93 

45.44 

43.1486 

7.14315 

29.1483 

57.1489 

2.2914 

47 

17JUL93 

69.08 

63.2411 

7.14315 

49.2408 

77.2414 

5.8389 

48 

18JUL93 

76.84 

74.5920 

7.14315 

60.5916 

88.5923 

2.2480 

49 

19JUL93 

52.32 

66. 6568 

7.14315 

52.6565 

80.6571 

-14.3368 

50 

20JUL93 

69.60 

58.8327 

7.14315 

44.8324 

72.8330 

10.7673 

51 

21JUL93 

88.04 

80.9423 

7.14315 

66.9420 

94.9426 

7.0977 

52 

22JUL93 

95.12 

91.9776 

7.14315 

77.9773 

105.978 

3.1424 

53 

23JUL93 

91.60 

94.9329 

7.14315 

80.9326 

108.933 

-3.3329 

54 

24JUL93 

81.40 

80.9557 

7.14315 

66.9554 

94.956 

0.4443 

55 

25JUL93 

57.20 

65.1737 

7.14315 

51.1734 

79.174 

-7.9737 

56 

26JUL93 

74.96 

77.9551 

7.14315 

63.9548 

91.955 

-2.9951 

57 

27JUL93 

74.67 

75.2995 

7.14315 

61.2992 

89.300 

-0.6295 

58 

28JUL93 

89.72 

80.4016 

7.14315 

66.4013 

94.402 

9.3184 

59 

29JUL93 

93.48 

79.8823 

7.14315 

65.8820 

93.883 

13.5977 

60 

30JUL93 

64.12 

64.5478 

7.14315 

50.5475 

78.548 

-0.4278 
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61 

31JUL93 

71.16 

71.4109 

62 

01AUG93 

82.04 

69.5280 

63 

02AUG93 

68.44 

68.1764 

64 

03AUG93 

48.87 

39.3754 

65 

04AUG93 

50.40 

47.8113 

66 

05AUG93 

55.88 

58.3703 

67 

06AUG93 

39.00 

31.7075 

68 

07AUG93 

35.84 

43.1916 

69 

08AUG93 

51.12 

48.4239 

70 

09AUG93 

69.64 

68.7672 

71 

10AUG93 

71.88 

71.0996 

72 

11AUG93 

76.55 

89.6603 

73 

12AUG93 

64.08 

72.6781 

74 

13AUG93 

52.08 

52.7281 

75 

14AUG93 

55.24 

49.7652 

76 

15AUG93 

63.16 

63.1766 

77 

16AUG93 

75.00 

69.7847 

78 

17AUG93 

77.52 

75.1357 

79 

18AUG93 

76.24 

65.5987 

80 

19AUG93 

101.64 

90.7582 

81 

20AUG93 

75.96 

76.8325 

82 

21AUG93 

74.12 

60.1109 

83 

22AUG93 

70.60 

68.3250 

84 

23AUG93 

93.84 

98.5485 

85 

24AUG93 

57.42 

53.9654 

86 

25AUG93 

77.28 

76.8137 

87 

26AUG93 

101.68 

94.4690 

88 

27AUG93 

67.52 

62.5746 

89 

28AUG93 

73.64 

74.3456 

90 

29AUG93 

70.88 

70.7708 

91 

30AUG93 

60.84 

76.0581 

92 

31AUG93 

75.96 

69.8248 

7.14315 

57.4106 

85.411 

-0 

7.14315 

55.5277 

83.528 

12 

7.14315 

54.1761 

82.177 

0 

7.14315 

25.3751 

53.376 

9 

7.14315 

33.8110 

61.812 

2 

7.14315 

44.3700 

72.371 

-2 

7.14315 

17.7072 

45.708 

7 

7.14315 

29.1913 

57.192 

-7 

7.14315 

34.4236 

62.424 

2 

7.14315 

54.7669 

82.767 

0 

7.14315 

57.0993 

85.100 

0 

7.14315 

75.6599 

103.661 

-13 

7.14315 

58.6777 

86.678 

-8, 

7.14315 

38.7277 

66.728 

-0, 

7.14315 

35.7649 

63.766 

5, 

7.14315 

49.1763 

77.177 

-0. 

7.14315 

55.7844 

83.785 

5, 

7.14315 

61.1354 

89.136 

2, 

7.14315 

51.5984 

79.599 

10. 

7.14315 

76.7579 

104.759 

10. 

7.14315 

62.8322 

90.833 

-0. 

7.14315 

46.1106 

74.111 

14. 

7.14315 

54.3247 

82.325 

2. 

7.14315 

84.5482 

112.549 

-4. 

7.14315 

39.9650 

67.966 

3. 

7.14315 

62.8134 

90.814 

0. 

7.14315 

80.4686 

108.469 

7. 

7.14315 

48.5743 

76.575 

4. 

7.14315 

60.3453 

88.346 

-0. 

7.14315 

56.7705 

84.771 

0. 

7.14315 

62.0578 

90.058 

-15. 

7.14315 

55.8245 

83.825 

6. 

.2509 

.5120 

.2636 

.4946 

.5887 

.4903 

.2925 

.3516 

.6961 

.8728 

.7804 

.1103 

.5981 

.6481 

.4748 

.0166 

.2153 

,3843 

.6413 

8818 

8725 

0091 

2750 

7085 

4546 

4663 

2110 

9454 

7056 

1092 

2181 

1352 
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Appendix  D.  Yorkville  1992  and  1993  Data 


/.  Yorkville  1992  meteorological  and  gas  data 


Date 

TEMP 

RELHUM 

pBARO 

SOLRAD 

OZONE 

inax03 

CO 

S02 

NO 

NOy 

6-1- 

-92 

15 

.30 

74. 

,80 

728 

.14 

408 

52 

46. 

46 

58 

262. 

,40 

14 

.14 

2. 

72 

16. 

98 

6-2- 

-92 

21 

.56 

58. 

.18 

727 

.31 

621 

08 

81. 

46 

94 

276. 

16 

2 

.76 

0. 

43 

14. 

49 

6-3- 

-92 

16 

.74 

98. 

.08 

724 

.26 

129 

84 

35. 

76 

39 

242. 

,44 

5 

.78 

0. 

58 

5. 

71 

6-4- 

-92 

20 

.91 

73, 

.02 

720 

.38 

434 

04 

57. 

04 

64 

194. 

20 

1 

.78 

0. 

24 

6. 

42 

6-5- 

-92 

22 

.68 

69. 

.72 

722 

.61 

664 

72 

55. 

34 

65 

206. 

33 

1 

.62 

0. 

61 

9. 

08 

6-6- 

-92 

26 

.29 

56. 

,03 

726 

.03 

785 

52 

65. 

40 

68 

241. 

00 

2 

.34 

0. 

30 

11. 

57 

6-7- 

-92 

26 

.60 

72. 

,08 

728 

.58 

668 

64 

49. 

70 

56 

239. 

00 

0 

.78 

0. 

27 

7. 

08 

6-8- 

-92 

19 

.60 

93. 

,10 

728 

.58 

376 

20 

18. 

81 

25 

211. 

36 

19 

.99 

3. 

42 

10. 

45 

6-9- 

-92 

20 

.62 

91. 

,25 

725 

.26 

501 

12 

24. 

18 

37 

0 

.56 

0. 

36 

1. 

86 

6-10- 

-92 

25. 

.90 

77. 

41 

725, 

.97 

630 

84 

48. 

88 

59 

122. 

63 

2 

.18 

0. 

37 

8. 

59 

6-11- 

■92 

24 

.20 

88. 

,08 

726, 

.41 

356 

64 

57. 

26 

72 

281. 

08 

4 

.04 

0. 

96 

13. 

35 

6-12- 

•92 

16, 

.56 

93. 

03 

728, 

.02 

150 

96 

44. 

92 

48 

293. 

72 

0 

.  19 

1. 

32 

15. 

47 

6-13- 

■92 

18, 

.30 

99. 

,10 

726, 

.08 

138 

88 

24. 

78 

33 

304. 

00 

0 

.46 

2. 

45 

15. 

00 

6-14- 

■92 

23, 

.63 

85. 

80 

725, 

.46 

472 

04 

40. 

47 

48 

181. 

16 

0 

.26 

0. 

21 

8. 

47 

6-15- 

■92 

26. 

.90 

73. 

78 

726. 

.97 

547 

38 

52. 

13 

58 

213. 

63 

0 

.57 

0. 

27 

6. 

04 

6-16- 

■92 

28. 

.49 

66. 

84 

729. 

.59 

591 

92 

81. 

83 

101 

224. 

88 

6 

.24 

0. 

33 

13. 

92 

6-17- 

■92 

26. 

.16 

78. 

13 

731. 

.48 

591 

44 

56. 

92 

63 

233. 

40 

0 

.41 

0. 

83 

12. 

46 

6-18- 

•92 

26. 

.58 

78. 

80 

729. 

.53 

471 

00 

41. 

67 

53 

210. 

00 

23 

.26 

3. 

10 

17. 

45 

6-19- 

■92 

27. 

.68 

65. 

77 

725. 

.23 

801 

64 

65. 

62 

77 

223. 

68 

4 

.17 

0. 

40 

11. 

96 

6-20- 

■92 

26. 

.51 

56. 

34 

723. 

.95 

792 

96 

67. 

74 

74 

203. 

40 

3 

.07 

0. 

32 

7. 

09 

6-21- 

■92 

23. 

.25 

73. 

94 

725. 

.81 

649 

84 

52. 

41 

64 

221. 

48 

1 

.31 

0. 

34 

5. 

55 

6-22- 

■92 

25. 

,08 

60. 

75 

727. 

,25 

687 

00 

68. 

88 

89 

193. 

50 

6 

.11 

2. 

20 

19. 

33 

6-23- 

■92 

25. 

.57 

55. 

41 

725. 

.66 

788 

60 

81. 

63 

107 

229. 

60 

8 

.35 

0. 

39 

12. 

02 

6-24- 

92 

22. 

,97 

48. 

76 

723. 

,77 

826 

48 

64. 

81 

71 

221. 

36 

1 

.23 

0. 

26 

8. 

02 

6-25- 

92 

23. 

,49 

58. 

78 

724, 

.32 

727 

68 

64. 

44 

70 

231. 

56 

2 

.40 

0. 

34 

8. 

19 

6-26- 

■92 

22. 

,18 

75. 

01 

724, 

,31 

600 

24 

58. 

90 

67 

245. 

24 

2 

.77 

0. 

40 

7. 

99 

6-27- 

92 

20. 

,84 

74. 

64 

724. 

.60 

638 

08 

57. 

88 

72 

257. 

88 

2 

.22 

0. 

50 

8. 

10 

6-28- 

■92 

19. 

,97 

75. 

82 

725, 

,03 

552 

36 

75. 

90 

87 

234. 

20 

5 

.67 

0. 

69 

14. 

11 

6-29- 

92 

26. 

,70 

61. 

76 

726. 

.10 

694 

60 

76. 

08 

83 

178. 

44 

10 

.13 

0. 

46 

14. 

23 

6-30- 

92 

23. 

55 

81. 

65 

725. 

,46 

4  63 

04 

46. 

46 

55 

126. 

04 

12 

.49 

0. 

99 

7. 

88 

7-1- 

92 

24. 

,12 

89. 

38 

726. 

.23 

504 

38 

33. 

26 

48 

162. 

21 

0 

.66 

0. 

42 

3. 

02 

7-2- 

92 

26. 

27 

80. 

74 

727. 

,55 

603 

32 

, 

. 

175. 

56 

0 

.24 

0. 

39 

4. 

86 

7-3- 

92 

24. 

00 

86. 

80 

729. 

,34 

37. 

81 

44 

204  . 

56 

0 

.21 

0. 

91 

5. 

36 

7-4- 

92 

27. 

65 

68. 

60 

731. 

.11 

56. 

01 

66 

187. 

44 

2 

.20 

0. 

30 

6. 

80 

7-5- 

92 

29. 

09 

62. 

20 

728. 

03 

41. 

43 

50 

146. 

44 

0 

.44 

0. 

23 

6. 

82 

7-6- 

92 

27. 

43 

77. 

48 

727. 

.36 

41. 

42 

47 

166. 

52 

1 

.20 

0. 

69 

5. 

19 

7-7- 

92 

30. 

38 

66. 

93 

730. 

60 

64. 

87 

83 

181. 

84 

8 

.55 

0. 

90 

12. 

47 

7-8- 

92 

31. 

86 

47. 

74 

732. 

.03 

. 

. 

167. 

92 

1 

.01 

1. 

67 

9. 

86 

7-9- 

92 

30. 

83 

61. 

07 

731. 

44 

48. 

70 

55 

189. 

28 

1 

.33 

0. 

38 

6. 

80 

7-10- 

92 

31. 

20 

64. 

24 

731. 

,12 

57. 

44 

64 

222. 

40 

2 

.71 

0. 

46 

8. 

69 

7-11- 

92 

31. 

28 

57. 

62 

730. 

83 

49. 

69 

59 

160. 

16 

1 

.48 

1. 

03 

9. 

32 

7-12- 

92 

33. 

03 

50. 

28 

730. 

74 

53. 

78 

67 

194. 

56 

4 

.44 

0. 

60 

9. 

12 

7-13- 

92 

30. 

06 

72. 

13 

730. 

66 

53. 

47 

61 

227. 

52 

4 

.14 

0. 

39 

8. 

38 

7-14- 

92 

29. 

39 

65. 

76 

729. 

68 

40. 

68 

50 

184. 

00 

2 

.46 

0. 

37 

5. 

78 

88 


7-15-92 

24.86 

88.26 

728.10 

. 

34.75 

45 

248.32 

2.76 

0.74 

5.06 

7-16-92 

27.28 

76.97 

729.67 

. 

35.67 

40 

217.12 

0.94 

0.35 

4.46 

7-17-92 

27.35 

77.54 

731.04 

. 

31.15 

39 

223.72 

2.38 

0.56 

4.79 

7-18-92 

24.22 

83.84 

731.72 

. 

39.78 

50 

229.64 

1.65 

0.42 

3.52 

7-19-92 

• 

64.78 

731.38 

. 

56.23 

72 

268.48 

27.36 

10.04 

37.54 

7-20-92 

. 

60.31 

. 

. 

74.16 

81 

301.24 

3.78 

0.27 

13.76 

7-21-92 

25.96 

66.75 

730.50 

. 

49.04 

54 

211.54 

3.30 

0.33 

9.20 

7-22-92 

23.46 

90.78 

730.38 

. 

35.71 

41 

205.09 

0.45 

1.17 

5.28 

7-23-92 

26.22 

80.08 

731.90 

. 

31.52 

37 

170.72 

1.01 

0.45 

4.19 

7-24-92 

26.97 

80.40 

733.37 

. 

33.04 

41 

161.28 

2.08 

0.49 

4.93 

7-25-92 

22.28 

67.26 

732.71 

. 

40.65 

50 

151.68 

2.72 

0.65 

7.97 

7-26-92 

26.57 

71.92 

729.36 

. 

38.93 

50 

166.84 

3.04 

0.72 

8.78 

7-27-92 

25.02 

81.18 

727.30 

. 

46.36 

54 

182.16 

0.78 

0.39 

5.85 

7-28-92 

24.80 

72.08 

728.75 

542.92 

48.56 

62 

208.72 

0.94 

0.55 

7.14 

7-29-92 

27.27 

67.07 

728.72 

655.56 

57.61 

67 

189.24 

1.20 

0.38 

8.20 

7-30-92 

26.05 

80.49 

729.14 

613.96 

47.28 

54 

160.44 

1.76 

0.51 

7.34 

7-31-92 

26.09 

79.96 

729.03 

507.84 

37.46 

43 

130.44 

1.40 

0.40 

4.61 

8-1-92 

23.82 

69.86 

731.04 

727.80 

51.31 

57 

171.76 

1.98 

0.44 

7.94 

8-2-92 

25.05 

57.91 

730.36 

772.20 

52.54 

55 

153.76 

0.98 

0.30 

7.41 

8-3-92 

22.28 

56.95 

728.65 

755.48 

60.02 

68 

165.16 

7.15 

0.64 

8.92 

8-4-92 

22.40 

69.94 

728.38 

617.36 

53.06 

57 

176.72 

0.80 

0.31 

7.24 

8-5-92 

26.10 

62.88 

730.31 

748.48 

61.82 

73 

195.72 

2.42 

0.35 

9.42 

8-6-92 

24.10 

80.56 

731.25 

382.64 

85.78 

113 

344.17 

8.40 

4.71 

42.33 

8-7-92 

23.23 

90.05 

731.96 

376.92 

44.76 

60 

227.32 

0.97 

0.61 

7.08 

8-8-92 

26.52 

76.74 

732.60 

442.80 

55.21 

59 

152.92 

16.00 

0.30 

7.42 

8-9-92 

27.80 

71.90 

731.71 

619.48 

49.92 

57 

144.24 

4.61 

0.17 

5.43 

8-10-92 

29.11 

67.52 

729.12 

748.91 

54.35 

60 

201.43 

2.20 

0.19 

6.27 

8-11-92 

27.67 

62.06 

728.56 

752.48 

48.51 

59 

274.92 

3.18 

0.22 

5.10 

8-12-92 

25.73 

83.02 

729.32 

321.60 

48.66 

69 

. 

11.50 

0.25 

5.72 

8-13-92 

18.68 

99.97 

728.44 

100.80 

35.30 

40 

160.08 

4.08 

0.22 

3.29 

8-14-92 

23.54 

72.92 

729.36 

657.92 

39.30 

59 

134.23 

35.24 

4.50 

13.84 

8-15-92 

22.52 

78.65 

729.25 

346.33 

56.81 

72 

178.24 

15.77 

0.39 

7.27 

8-16-92 

21.35 

90.54 

729.62 

318.56 

48.20 

62 

291.84 

1.64 

0.32 

5.50 

8-17-92 

. 

86.49 

729.94 

442.72 

39.46 

55 

211.12 

1.38 

0.47 

4.83 

8-18-92 

. 

77.69 

731.13 

565.60 

51.76 

66 

204.92 

3.24 

0.28 

4.77 

8-19-92 

24.39 

65.05 

729.62 

476.04 

69.52 

74 

222.00 

9.82 

0.18 

6.76 

8-20-92 

24.55 

63.26 

728.85 

560.92 

73.20 

91 

267.20 

9.18 

0.18 

8.97 

8-21-92 

. 

100.00 

731.18 

83.60 

25.76 

39 

353.12 

4.38 

1.21 

11.62 

8-22-92 

23.48 

88.75 

731.94 

288.00 

33.71 

42 

206.80 

14.50 

0.52 

6.78 

8-23-92 

24.34 

86.04 

733.37 

477.08 

34.42 

41 

211.88 

4.96 

0.45 

6.07 

8-24-92 

25.18 

77.07 

734.79 

577.60 

55.33 

66 

155.16 

2.18 

0.68 

8.32 

8-25-92 

21.30 

71.89 

734.08 

530.32 

62.43 

76 

170.84 

6.12 

0.75 

9.32 

8-26-92 

26.59 

81.26 

732.20 

463.48 

26.35 

36 

189.35 

9.70 

0.58 

5.62 

8-27-92 

23.10 

90.58 

725.81 

301.64 

23.83 

27 

160.28 

6.80 

0.76 

4.11 

8-28-92 

20.02 

85.44 

727.03 

444.40 

25.81 

41 

164.57 

2.62 

0.65 

3.74 

8-29-92 

20.60 

68.97 

730.14 

759.36 

41.32 

55 

139.00 

22.09 

3.65 

12.03 

8-30-92 

23.78 

64.68 

731.82 

634.16 

57.32 

73 

135.44 

3.39 

0.14 

4.83 

8-31-92 

24.51 

64.29 

732.08 

674.20 

76.30 

86 

131.72 

0.31 

10.84 
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II.  Yorkville  1993  meteorological  and  gas  data 


York 

TEMP 

RELHUM 

pBARO 

SOLRAD 

OZONE 

max03 

CO 

S02 

NO 

NOy 

6-1-93 

. 

. 

. 

902.13 

55.04, 

63 

293.16 

2.08 

0.38 

4.07 

6-2-93 

. 

. 

. 

794.88 

68.24 

83 

374.60 

4.96 

0.60 

6.99 

6-3-93 

26.59 

56.94 

723.66 

869.68 

71.48 

83 

300.96 

3.14 

0.17 

4.44 

6-4-93 

26.54 

72.32 

723.67 

786.00 

36.20 

39 

277.80 

0.62 

0.27 

2.82 

6-5-93 

24.85 

87.38 

726.75 

397.38 

44.04 

50 

315.67 

0.80 

0.24 

2.96 

6-6-93 

23.34 

56.66 

730.05 

877.32 

57.28 

73 

289.44 

2.38 

0.22 

3.81 

6-7-93 

28.64 

64.48 

729.54 

793.72 

62.52 

72 

344.84 

6.12 

0.26 

4.49 

6-8-93 

29.27 

63.20 

727.98 

806.48 

50.84 

57 

280.20 

2.45 

0.29 

3.41 

6-9-93 

29.00 

61.16 

728.86 

828.88 

42.20 

54 

246.12 

1.73 

0.25 

2.82 

6-10-93 

28.63 

58.81 

729.78 

790.00 

53.80 

61 

271.32 

2.71 

0.23 

3.27 

6-11-93 

28.72 

58.33 

729.63 

844.44 

52.28 

55 

280.68 

1.31 

0.22 

3.08 

6-12-93 

27.69 

64.90 

729.14 

746.16 

84.46 

119 

344.60 

32.20 

0.72 

11.04 

6-13-93 

25.91 

72.63 

729.07 

707.04 

71.88 

75 

373.00 

4.62 

0.24 

5.64 

6-14-93 

21.48 

91.27 

729.12 

371.16 

57.32 

74 

502.32 

2.14 

0.73 

7.49 

6-15-93 

25.39 

71.23 

729.31 

583.50 

65.12 

70 

355.56 

1.50 

0.24 

4.78 

6-16-93 

27.44 

64.90 

730.18 

738.88 

81.56 

97 

351.28 

13.89 

0.34 

8.07 

6-17-93 

27.38 

59.50 

732.74 

840.12 

76.84 

110 

337.80 

8.09 

0.48 

6.63 

6-18-93 

27.22 

63.35 

733.10 

753.40 

77.20 

100 

346.76 

2.51 

0.50 

7.04 

6-19-93 

26.41 

59.20 

731.54 

749.68 

73.36 

92 

305.44 

3.18 

0.25 

6.82 

6-20-93 

26.79 

56.45 

730.73 

590.76 

58.52 

70 

255.16 

22.94 

0.50 

6.80 

6-21-93 

20.04 

97.39 

729.50 

284.00 

36.36 

42 

315.12 

3.30 

0.47 

4.75 

6-22-93 

23.70 

82.71 

726.77 

429.43 

32.67 

42 

280.16 

1.13 

0.27 

2.78 

6-23-93 

27.82 

61.50 

727.56 

811.44 

93.76 

104 

423.20 

5.03 

0.32 

9.53 

6-24-93 

26.24 

74.95 

730.07 

560.64 

65.64 

80 

398.80 

1.25 

0.41 

7.04 

6-25-93 

23.11 

87.56 

731.17 

. 

41.40 

57 

327.00 

23.64 

1.29 

7.71 

6-26-93 

25.69 

76.17 

729.15 

. 

44.52 

56 

300.40 

0.94 

0.27 

3.18 

6-27-93 

27.35 

67.64 
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